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INTRODUCTION 

Epidemiologic  studies  have  highlighted  the  relationship  between  hormones  and  carcinogenesis  in 
the  breast  by  identifying  endocrine  risk  factors  for  breast  cancer  that  are  related  to  the  timing  of 
reproductive  endocrine  events,  such  as  menarche  and  menopause.  This  relationship  is  illustrated  by  the 
observation  that  women  who  undergo  an  early  first  full-term  pregnancy  have  a  significantly  reduced 
lifetime  risk  of  breast  cancer.  The  recognition  that  specific  reproductive  endocrine  events  alter  breast 
cancer  risk  in  a  predictable  fashion  raises  the  possibility  that  events  associated  with  a  decrease  in  breast 
cancer  risk,  such  as  early  first  full-term  pregnancy,  might  be  mimicked  pharmacologically.  As  such, 
understanding  the  mechanisms  by  which  these  events  influence  breast  cancer  risk  would  facilitate  the 
design  of  safe  and  effective  hormonal  chemoprevention  regimens.  In  addition,  the  testing  of  such 
regimens  would  be  facilitated  by  the  identification  of  biomarkers  that  accurately  reflect  early  biological 
changes  in  the  breast  associated  with  reproductive  endocrine  events  that  alter  breast  caneer 
susceptibility.  Since  the  impact  of  reproductive  endocrine  history  both  on  mammary  gland  development 
and  on  breast  cancer  risk  most  likely  results  from  differential  effects  of  hormones  on  particular 
subpopulations  of  cells,  understanding  these  changes  will  almost  certainly  require  a  thorough 
understanding  of  the  cell  types  present  in  the  breast  and  the  manner  in  which  hormones  affect  their 
normal  programs  of  differentiation  and  development.  Unfortunately,  pursuit  of  this  goal  has  been 
hampered  by  the  fact  that  few,  if  any,  molecular  biomarkers  have  been  identified  to  date  that  are  specific 
for  changes  in  the  breast  that  occur  as  a  result  of  reproductive  endocrine  events  known  to  influence 
breast  cancer  risk.  Such  biomarkers  are  essential  both  for  understanding  the  impact  of  reproductive 
history  on  mammary  epithelial  biology  in  general,  and  the  molecular  and  cellular  basis  of  parity-induced 
protection  against  breast  cancer  in  particular,  as  well  as  for  the  rational  design  and  testing  of  hormonal 
chemoprevention  regimens  aimed  at  mimicking  this  naturally  occurring  protective  event. 

We  hypothesize  that  an  early  first  full-term  pregnancy  results  in  a  permanent  change  in  the  breast 
that  confers  a  decreased  risk  for  the  subsequent  development  of  breast  cancer.  We  further  hypothesize 
that  parity-induced  changes  in  the  number  and  distribution  of  specific  epithelial  cell  subtypes  in  the 
breast  may  contribute  to  parity-induced  changes  in  breast  cancer  risk.  To  address  these  hypotheses,  the 
specific  aims  of  this  proposal  were  designed  to  identify  and  evaluate  molecular  biomarkers  for  parity- 
induced  changes  in  the  breast  by  isolating  genes  that  are  differentially  expressed  between  the  parous  and 
nulliparous  rodent  breast  using  DNA  chip  technology.  Each  of  the  technical  objectives  and  tasks 
proposed  have  been  completed  on  schedule. 

In  this  report  we  have  used  high-density  oligonucleotide  microarrays  to  analyze  the  impact  of 
early  first  full-term  pregnancy  on  global  gene  expression  profiles  within  the  murine  mammary  gland. 
This  approach  has  led  to  the  identification  of  a  panel  of  genes  whose  expression  in  the  mammary  gland 
is  persistently  altered  as  a  consequence  of  parity  in  multiple  strains  of  mice  as  well  as  in  two  widely 
used  rat  models  for  parity-induced  protection  against  breast  cancer.  Our  findings  demonstrate  that  parity 
induces  the  persistent  downregulation  of  multiple  genes  encoding  epithelial  growth  factors  as  well  as  the 
persistent  upregulation  of  the  growth-inhibitory  molecule,  TGF-J33,  and  several  of  its  downstream 
targets.  Our  findings  further  indicate  that  parity  results  in  a  persistent  increase  in  the  differentiated  state 
of  the  mammary  gland  as  well  as  permanent  changes  in  the  hematopoietic  cell  types  resident  within  the 
gland.  These  findings  define  a  developmental  state  of  the  mammary  gland  that  is  refractory  to 
carcinogenesis  and  suggest  novel  hypotheses  for  mechanisms  by  which  parity  may  modulate  breast 
cancer  risk.  We  believe  that  these  studies  represent  an  important  step  towards  determining  the 
mechanisms  by  which  breast  cancer  susceptibility  is  modulated  by  reproductive  history. 
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BODY 

At  the  outset  of  this  study,  we  hypothesized  that  generating  a  sufficiently  robust  panel  of  parous-  and 
nulliparous-specific  markers  would  require  a  powerful  parallel  screening  technique  that  would  permit 
the  sensitive,  reliable  and  reproducible  detection  of  differentially  expressed  genes.  To  accomplish  this 
aim,  we  chose  to  determine  gene  expression  levels  by  hybridization  of  mRNA  to  high-density 
oligonucleotide  arrays.  This  approach  has  turned  out  to  be  extremely  successful.  While  several  of  the 
highlights  of  our  experimental  studies  are  presented  in  this  report,  the  complete  results  obtained  are 
described  in  the  accompanying  manuscript  from  my  laboratory  entitled,“Early  first  full-term  pregnancy 
results  in  permanent  changes  in  gene  expression  in  the  murine  mammary  gland”  by  D’Cruz  et  al.  This 
manuscript  has  been  submitted  for  publication.  Please  refer  to  this  manuscript  for  all  details. 

Technical  Objective  I:  Isolate  genes  expressed  in  the  breast  in  a  parity-specific  manner. 

To  determine  whether  our  model  system  behaved  similarly  to  other  established  systems  for  parity- 
induced  protection  against  breast  cancer,  we  initially  compared  the  morphology  of  mammary  glands 
from  rodents  that  have  undergone  a  single-round  of  pregnancy,  21  days  of  lactation  and  4  weeks  of 
regression  to  the  morphology  of  mammary  glands  from  nulliparous  littermates  (Fig.  1).  It  has 
previously  been  noted  that  mammary  glands  of  parous  woman  have  more  differentiated  lobules  than 
seen  in  nulliparous  woman[l].  Whole  mount  analysis  of  mammary  glands  from  virgin  and  parous 


Rat  C57B1/6  FVB  FVB-  H&E 


Virgin 


Parous 


Figure  1:  Parity-induced  changes  in  mammary  gland  morphology 


Sprague-Dawley  rats  demonstrate  marked  differences  in  the  architecture  of  the  virgin  and  parous  glands. 
The  parous  gland  is  more  highly  branched,  suggesting  that  it  is  permanently  altered  by  parity-related 
changes.  We  have  also  examined  two  independent  strains  of  mice,  C57BL/6  and  FVB  to  determine 
whether  these  morphological  features  appear  in  the  murine  system  as  well.  Parous  animals  were  mated 
at  4  wks  of  age,  a  time  corresponding  to  puberty  in  order  to  approximate  “early-parity”.  Although  the 
virgin  murine  glands  appear  less  developed  as  compared  to  the  rat,  similar  overall  differences  are 
observed  between  the  virgin  and  parous  glands  of  mice  and  rats.  Again  the  parous  gland  is  more  highly 
branched  (Fig  1).  Analysis  of  H&F  stained  glands  demonstrate  that  there  appears  to  be  roughly  the 
same  amount  of  epithelium,  an  observation  that  has  been  confirmed  by  quantitation  of  epithelium  by 
Cytokeratin  18  expression  levels.  The  structural  changes  which  accompany  parity  in  the  mammary 
gland  appear  to  be  permanent  as  the  differences  in  the  gland  remain  for  essentially  the  life  of  the  animal 
(data  not  shown). 
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Changes  in  the  breast  that  are  induced  by  reproductive  events  such  as  parity  are  likely  to  involve 
complex  changes  in  the  expression  of  multiple  genes.  Many  techniques  for  comparing  mRNA 
populations,  such  as  differential  display,  subtractive  hybridization,  and  serial  analysis  of  gene  expression 
(SAGE)  have  proven  to  be  labor-intensive  and  relatively  insensitive,  particularly  for  genes  expressed  at 
very  low  levels.  In  addition,  these  techniques  typically  generate  high  percentages  of  false-positive 
clones.  As  such,  we  have  chosen  a  technical  approach  that  permits  the  direct  assessment  of  mRNA 
expression  profiles  from  large  numbers  of  genes  in  parallel.  This  approach,  referred  to  as  Affymetrix 
GeneChip  technology,  involves  the  determination  of  the  relative  abundance  of  mRNAs  based  on 
hybridization  of  complex  mRNA  populations  to  high-density  arrays,  or  "chips"  containing  arrayed 
synthetic  oligonucleotides [2-4].  Each  chip  is  composed  of  65,000  different  oligonucleotides  of  defined 
sequence  which  are  synthesized  directly  on  derivatized  glass  slides.  Each  gene  is  represented  on  the 
array  by  20  specific,  unique  25-mer  oligonucleotide  probe  pairs,  such  that  each  chip  permits  the 
detection  of  approximately  1,600  genes.  A  combination  of  arrays  permits  expression  levels  to  be 
determined  for  more  than  6,800  full-length  genes  in  the  UniGene,  GenBank  and  TIGR  databases. 

mRNA  isolated  from  a  sample  of  interest  is  used  to  generate  an  antisense  poly(A)"''  probe  with 
incorporated  biotin-labeled  nucleotides.  This  mRNA  probe  is  then  hybridized  to  the  array  using  the 
GeneChip  fluidics  station  which  automates  sample  introduction,  hybridization,  incubation,  and  washing. 
Quantitative  detection  of  mRNA  bound  to  each  oligonucleotide  cell  is  achieved  by  measuring 
fluorescence  of  bound  phycoerythrin-streptavidin  using  a  Hewlett-Packard  GeneArray  scanner.  Finally, 
the  GeneChip  operating  system  software  automatically  calculates  intensity  values  for  each  probe  cell. 
Expression  for  each  gene  is  determined  by  analyzing  the  extent  of  binding  to  the  20  pairs  of 
oligonucleotides  specific  for  that  gene.  Fluorescence  intensity  from  the  hybridized  arrays  directly 
correlate  with  mRNA  abundance  levels.  The  GeneChip  instrumentation  system  necessary  for  these 
experiments  has  recently  been  purchased  from  Affymetrix,  Inc.  by  the  Institute  for  Human  Gene 
Therapy.  This  equipment  will  be  located  within  the  Chodosh  laboratory. 

GeneChip  technology  permits  the  quantitative,  sensitive,  specific  and  reproducible  determination  of 
gene  expression  profiles  [2-4].  A  high  percentage  of  clones  identified  are  true  positives,  thereby 
minimizing  subsequent  labor  invested  in  pursuing  false  leads.  Differences  in  expression  levels  of 
approximately  2  to  2.5-fold  can  be  reliably  distinguished.  The  detection  limit  is  approximately  2-5 
copies  per  cell,  and  is  quantitative  over  three  orders  of  magnitude  [2-4].  Each  array  contains  probes 
comlementary  to  multiple  reference  genes  that  permit  data  within  and  among  experiments  to  be 
normalized  and  quantitated  [2-4].  In  addition,  multiple  internal  referenee  points  permit  the 
normalization  and  quantitation  of  data  from  different  experiments.  Repeat  hybridization  with  the  same 
probe  is  reported  to  yield  only  one  in  3,000  genes  that  differ  by  more  than  3-fold  [2-4].  Hybridization 
with  probes  independently  prepared  by  different  operators  from  the  same  sample  source  is  reported  to 
yield  less  than  1  in  1,000  genes  that  differ  by  more  than  a  factor  of  three.  Compared  to  array -based 
methods  that  involve  the  spotting  of  cDNAs  onto  nylon  membranes,  an  additional  advantage  of  the 
oligonucleotide-based  approach  is  that  highly  specific  binding  and  detection  can  be  achieved  even 
within  conserved  gene  families.  Finally,  this  parallel  approach  lends  itself  to  increases  in  scale  to  permit 
large  numbers  of  genes  to  be  monitored  and  is  thereby  ideally  suited  for  the  extensive  comparison  of 
complex  mRNA  populations  such  as  were  expected  in  this  proposal. 


Task  1:  Months  1-12:  Generation  of  mRNA  samples  representing  reproductive  endpoints. 

This  task  has  been  completed  on  schedule.  Inbred  FVB  mice  (which  do  not  harbor  active  MMTV 
pro  virus)  were  bred  and  maintained  in  a  barrier  facility  with  a  regular  12  hour  light/12  hour  dark  cycle  to 
facilitate  regular  estrous  cycling.  Animals  were  given  food  and  water  ad  libitum.  Mice  were  housed  4 
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per  cage  during  virgin  development  and  postlactational  regression,  and  individually  during  pregnancy 
and  lactation.  Nulliparous  controls  were  housed  identically  to  parous  animals.  Daily  vaginal  smears 
were  performed  on  animals  to  document  regular  cycling.  Total  RNA  was  isolated  from  the  mammary 
glands  of  mice  that  had  either  undergone  an  early  first  full-term  pregnancy  followed  by  3  weeks  of 
lactation  and  4  weeks  of  postlactational  regression,  or  that  were  maintained  as  age-matched  nulliparous 
controls.  Animals  were  euthanized  by  carbon  dioxide  asphyxiation  at  specified  developmental  time 
points,  and  the  number  4,  5  and  6  mammary  gland  pairs  were  harvested  and  snap  frozen  in  liquid 
nitrogen,  following  removal  of  the  lymph  node  embedded  in  gland  4.  Mammary  tissue  from  each 
animal  was  also  fixed  in  10%  neutral  buffered  formalin  and  embedded  in  paraffin  for  in  situ 
hybridization  studies.  Total  RNA  was  isolated  from  tissues  by  homogenization  in  guanidinium 
isothiocyanate,  followed  by  ultracentrifugation  through  cesium  chloride.  In  total,  RNA  was  harvested 
from  parous  and  nulliparous  animals  in  three  independent  experiments  with  each  experiment  comprised 
of  20  parous  and  20  nulliparous  animals.  Additional  breeding  of  animals  is  ongoing  in  order  to  generate 
developmental  time  points  representing  other  reproductive  variables  such  as  early  vs.  late  first  full-term 
pregnancy  and  duration  of  postlactational  regression. 


Task  2:  Months  6-18:  GeneChip-based  screening. 

This  task  has  been  completed  on  schedule.  Standard  Affymetrix  protocols  were  used  to 
generate  and  hybridize  probes  prepared  from  nulliparous  and  parous  RNA  samples,  and  to  quantitate 
and  analyze  gene  expression  levels.  RNA  isolated  from  nulliparous  and  parous  animals  was  subject  to 
reverse  transcription  to  make  first  strand  cDNA  that  incorporated  a  T7  RNA  polymerase  promoter  at  the 
3’  end.  Following  completion  of  the  second  strand  of  cDNA,  this  material  was  used  in  an  in  vitro 
transcription  reaction  to  generate  biotin-labeled  cRNA  representing  each  experimental  pool  of  mammary 
gland  RNA.  This  in  vitro  transcribed  cRNA  was  then  fragmented  to  facilitate  hybridization. 
Hybridization  reactions  were  prepared  for  each  sample  using  fragmented  cRNA  in  addition  to 
standardized  internal  controls.  Hybridization  reactions  representing  each  of  the  3  nulliparous  samples  as 
well  as  each  of  the  3  parous  samples  were  hybridized  to  murine  6500  A,  B,  C  and  D  oligonuleotide 
arrays,  representing  approximately  5500  murine  genes. 

Using  GeneChip  software,  expression  levels  were  calculated  for  each  of  the  5500  genes 
represented  on  the  Mu6500  chips  for  all  six  samples.  This  information  was  imported  into  a  software 
database.  Confidence  intervals  were  generated  for  each  gene  in  order  to  identify  genes  that  were 
differentially  expressed  between  the  nulliparous  and  parous  states  at  a  confidence  level  >  90%.  In 
addition,  Affymetrix  algorithms  were  applied  to  each  individual  nulliparous/parous  pair  to  identify  those 
genes  that  whose  expression  was  considered  to  increase  or  decrease.  This  analysis  was  repeated  for 
each  of  the  three  experiments.  Genes  that  were  considered  to  either  increase  in  all  three  experiments  or 
decrease  in  all  three  experiments  were  identified.  This  approach,  combined  with  Northern  analysis  of 
candidate  genes  identified  by  microarray  analysis  as  being  differentially  expressed  in  two  out  of  three 
sets  of  nulliparous  and  parous  samples,  resulted  in  the  identification  of  41  differentially  expressed  genes 
(Table  1).  These  include  10  genes  that  are  preferentially  expressed  in  the  nulliparous  mammary  gland  as 
well  as  3 1  genes  that  are  preferentially  expressed  in  the  parous  mammary  gland.  Notably,  expression 
levels  of  the  epithelial  cell  marker,  cytokeratinlS,  do  not  differ  between  the  nulliparous  and  parous 
mammary  gland,  indicating  that  the  observed  changes  in  gene  expression  are  not  merely  a  consequence 
of  the  expansion  or  contraction  of  the  epithelial  cell  compartment  (Fig.  2). 
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Table  1;  Differentially  Expressed  Genes 


Genes  Down-Regulated 


Sequence  Identity 

Accession  Number 

Function 

Avg.  Fold  Change 

Northern 

Amphiregulin 

L41352 

Growth  Factor 

2.8 

Y 

Pleiotrophin 

D90225 

Growth  Factor 

2.0 

Y 

Insulin-Like  Growth  Factor  IB 

W 10072 

Growth  Factor 

2.2 

hD 

Insulin-Like  Growth  Factor  1A 

X04480 

Growth  Factor 

1.6 

Y 

Thyroid  Stimulating  Hormone  Receptor 

U02602 

Hormone  Receptor 

2.0 

Y 

Leptin 

U18812 

Hormone 

2.0 

Y 

Ig  Superfamily  Containing  Leucine-rich  repeat 

AA059664 

Cell  Adhesion 

2.5 

Y 

MUC  18 

AA088962 

Cell  Adhesion 

2.0 

Y 

Superoxide  Dismutase  111  (SODS) 

X84940 

Oxidoreductase 

2.5 

Y 

Carbonic  Anhydrase  III 

M27796 

Hydratase 

1.4 

Y 

Genes  Up-Regulated 


Sequence  Identity 

Accession  Number 

Function 

Avg.  Fold  Change 

Northern 

Whey  Addle  Protein 

J00544 

Dlfferentiation/Milk  Protein 

24  ' 

Y 

Gamma-Casein 

D10215 

Dlfferentiation/Milk  Protein 

20 

Y 

Aipha-Lactalbumin 

M80909 

Differentiation/Milk  Protein 

13 

Y 

Lactoferrin 

J03298 

DIfferentlation/lron  Transport 

4.8 

Y 

Alpha-Casein 

M36780 

Dlfferentiation/Milk  Protein 

4.2 

ISD 

WDNM1 

X93037 

Differentiation/Protease  Inhibitor 

3.5 

Y 

Beta-Casein 

X04490 

Differentiation/Milk  Protein 

3.3 

Y 

Kappa-Casein 

M10114 

Dlfferentiation/Milk  Protein 

2.9 

Y 

Adipocyte  Differentiation  Related  Protein 

M93275 

Differentiation 

1.8 

Y 

Carboxyl  Ester  Lipase 

U37386 

Lipid  Degradation 

6.6 

Y 

LPS-BInding  Protein 

X99347 

Antibacterial/Milk 

6.5 

Y 

Lysozyme  P 

M21050 

Antibacterial/Milk 

2.0 

kD 

Immunoalobulln  M  Heavy  Chain 

ET61785 

Immunoglobulin 

30 

Y 

Immunoglobulin  G  Heavy  Chain 

ET61798 . 

Immunoglobulin 

13 

Y 

Immunoglobulin  A  Heavy  Chain 

J00475 

Immunoglobulin 

21 

Y 

Kappa  Light  Chain 

X16678 

Immunoglobulin 

7.0 

Y 

Macrophage  Metalloelastase 

M82831 

Metalloprotease 

2.1 

Y 

Macrophage  Expressed  gene  1  (Mpegl) 

L20315 

Cell  Signalling 

1.7 

Y 

Lipocalin  2 

W13166 

Cell  Signalling 

2.6 

Y 

T-cell  death  associated  gene  (Tdag) 

U44088 

Cell  Signalling 

2.3 

Y 

Transforming  Growth  Factor  Beta-3 

M32745 

Growth  Inhibition 

1.9 

Y 

Clusterin 

L08325 

Apoptosis/TGF-B  Pathway 

2.0 

Y 

Eta-1 

X16151 

Cell  Signalling/TGF-B  Pathway 

8.8 

Y 

Id-2 

M69293 

Cell  Cycte/TGF-B  Pathway 

2.3 

Y 

Carbonic  Anhydrase  II 

K0081 1 

Hydratase 

3.0 

Y 

Cyclin  D1 

A1849928 

Cell  Cycle 

1.8 

Y 

CRBP1 

X60367 

Retinol-Binding 

3.5 

Y 

Connextn  26 

M81445 

Gap  Junction 

3.1 

Y 

Folate-Binding  Protein  1 

ET63126 

Receptor 

2.6 

Y 

Adenosine  Deaminase 

M10319 

Nucleotide  Metabolism 

3.0 

Y 

Chitinase  3-like  1  (BRP-39) 

X93035 

Glycoprotein 

3.2 

Y 

Table  1:  Genes  found  to  exhibit  a  persistent  change  in  expression  in  fuiiy  involuted  compared  to  age-matched 
nulliparous  mammary  glands.  Genes  whose  differential  expression  was  confirmed  by  Northern  analysis  are 
indicated  (Y).  Genes  whose  differential  expression  was  not  tested  (ND)  by  Northern  hybridization  were  listed 
if  they  displayed  a  consistent  change  in  regulation  in  3/3  microarray  experiments. 


Task  3:  Months  6-20:  Testing  of  candidate  genes  by  Northern  hybridization  and  RNase  protection. 

This  task  has  been  completed  on  schedule.  In  order  to  test  the  differential  expression  of  the 
candidate  genes  that  we  had  identified,  we  generated  400-600  base  pair  fragments  from  the  designated 
genes  using  RT-PCR  and  probed  Northern  blots  containing  one  virgin  and  one  parous  sample  identical 
to  material  that  was  used  in  the  GeneChip  analysis,  as  well  as  two  independent  virgin  and  parous _ 
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samples  to  control  for  biological  variability.  This  approach  confirmed  the  differential  expression  of  all 
of  the  genes  in  which  expression  either  increased  or  decreased  in  all  three  parous/nulliparous  sets,  as 
well  as  many  of  the  genes  whose  expression  either  increased  or  decreased  in  two  of  the  three 
parous/nulliparous  sample  sets,  that  also  had  a  confidence  level  >  90%  (Fig.  2).  Overall,  the  true¬ 
positive  rate  ascertained  to  date  is  approximately  70%.  This  true-positive  rate  is  considerably  higher 
than  that  from  any  of  the  procedures  that  we  have  previously  used  to  screen  for  differentially  expressed 
genes.  Approximately  80  genes  were  tested  by  Northern  hybridization. 
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Figure  2:  Differential  expression  of  genes  isolated  by  oligonucleotide  microarrays 


Since  the  phenomenon  of  parity-induced  protection  against  breast  cancer  is  conserved  among 
humans,  rats,  and  mice,  we  predicted  that  the  molecular  changes  that  underlie  this  effect  would  be 
conserved.  As  described  above,  microarray  analysis  of  gene  expression  profiles  in  the  mammary  glands 
of  parous  and  nulliparous  FVB  mice,  Balb/c  mice,  and  Lewis  rats  indicated  that  parity-dependent 
changes  in  gene  expression  are  largely  conserved  in  different  strains  of  mice  and  rats.  To  extend  these 
findings,  we  performed  Northern  analysis  on  pools  of  mammary  gland  total  RNA  isolated  from 
nulliparous  and  parous  129SvEv  and  Balb/c  mice,  as  well  as  from  nulliparous  and  parous  Sprague- 
Dawley  rats,  another  widely  used  model  for  parity-induced  protection  against  breast  cancer.  We  found 
that  all  ten  genes  examined  in  129SvEv  and  Balb/c  mouse  strains  were  differentially  expressed  in  a 
manner  similar  to  that  observed  in  FVB  mice  (Fig.  3B  and  data  not  shown).  Likewise,  each  of  the  nine 
genes  examined  in  the  rat  exhibited  a  parity-dependent  differential  pattern  of  expression  identical  to  that 
observed  in  the  mouse  (Fig.  3C  and  data  not  shown).  When  taken  together  with  the  above  microarray 
data,  these  findings  indicate  that  many  of  the  parity-related  molecular  changes  identified  in  this  study  are 
conserved  among  species  that  display  parity-induced  protection  against  mammary  tumorigenesis.  (Fig. 
3). 
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Figure  3:  Differential  expression  of  selected  genes  in  mice  and  rats 


Technical  Objective  II:  Analyze  genes  expressed  in  a  developmental  stage-specific  manner. 

Task  1:  Months  6-24:  Determine  the  developmental  expression  pattern  ofcDNA  markers. 

This  task  has  been  completed  on  schedule.  We  have  used  the  genes  that  we  have  identified  as 
being  differentially  expressed  to  define  parity-induced  cellular  and  molecular  changes  that  occur  in  the 
breastby  examining  the  developmental  expression  pattern  for  those  genes  that  we  have  confirmed  as 
being  differentially  expressed  as  a  function  of  parity.  Gene  expression  patterns  have  been  determined 
during  mammary  gland  development  using  Northern  hybridization.  Developmental  time  points 
analyzed  included  virgin  development  from  a  point  preceding  puberty  through  the  completion  of  ductal 
morphogenesis  (2 , 5  ,  10  ,  and  15  weeks);  early,  mid  and  late  pregnancy  (day  7,  14,  and  21);  mid 
lactation  (day  9),  and  early,  mid  and  late  postlactational  regression  (day  2,  7,  and  28).  For  complete 
evaluation,  see  accompanying  manuscript.  As  examples,  the  developmental  expression  of  the  epidermal 
growth  factor,  amphiregulin,  and  the  heparin-binding  mitogen,  pleiotrophin,  are  shown  in  Figure  4.  In 
adult  animals,  both  amphiregulin  and  pleiotrophin  are  expressed  in  a  nulliparous-specific  manner. 
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Developmentally,  expression  of  each  of  these  genes  appeared  to  be  induced  during  puberty  and  to 
remain  elevated  thereafter  in  the  virgin  gland.  Interestingly,  expression  of  both  genes  subsequently 
decreases  during  mid-pregnancy  and  remains 
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Fig.  4:  Down-regulation  of  amphiregulin  and  pleiotrophin  expression  as  a  function  of  parity 


suppressed  during  lactation  and  postlactational  regression. 

The  developmental  expression  patterns  of  numerous  other  genes  identified  in  this  screen  have 
also  been  examined.  For  example,  a  second  category  of  genes  whose  expression  was  elevated  in  the 
parous  involuted  mammary  gland  included  those  that  are  specifically  expressed  in  hematopoietic  cells 
such  as  B-lymphocytes  {kappa  light  chain  and  the  IgG,  IgM,  and  IgA  heavy  chains),  T-lymphocytes  {T- 
cell  death  associated  gene),  and  macrophages  {macrophage  metalloelastase  and  macrophage  expressed 
Lipocalin-2,  that  are  either  expressed  by  hematopoietic  cells  or  are  chemoattractants  for  these  cell  types 
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Fig.  5:  Up-regulation  of  hematopoietic  cell  gene  expression  as  a  function  of  parity 


gene  1).  Additional  genes,  including  Early  T-cell  activation  protein  {Eta-1),  EPS-Binding  protein  and 
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were  also  persistently  upregulated  as  a  consequence  of  parity.  The  marked  upregulation  of  these  genes 
in  the  parous  gland  suggests  that  cells  of  B-lymphocyte,  T-lymphocyte,  and  macrophage  lineages  may 
be  more  abundant  in  the  parous  mammary  gland.  To  investigate  this  hypothesis  further,  we  examined 
the  expression  patterns  for  several  immune-related  genes  during  stages  of  postnatal  mammary 
development  representing  puberty,  pregnancy,  lactation  and  involution.  Northern  analysis  revealed  that 
expression  of  the  B-cell-specific  gene,  K-light  chain  (k-LQ,  is  first  detected  during  lactation  with 
elevated  levels  of  expression  persisting  throughout  postlactational  involution  (Fig.  5).  Macrophage 
metalloelastase  (MME  or  MMP-12)  was  also  found  to  be  persistently  upregulated  in  the  mammary  gland 
as  a  consequence  of  parity  (Fig.  5A).  MME  is  a  secreted  metalloprotease  that  cleaves  plasminogen  to 
generate  angiostatin,  a  potent  inhibitor  of  endothelial  cell  proliferation.  Northern  hybridization 
demonstrated  that  MME  expression  levels  increase  dramatically  at  day  7  of  postlactational  involution 
and  remain  elevated  compared  to  age-matched  nulliparous  controls  following  28  days  of  involution  (Fig. 
5).  Finally,  we  examined  the  developmental  basis  for  the  persistent  parity-dependent  upregulation  of 
Eta-1  expression.  Northern  analysis  revealed  that  Eta-1  expression,  which  has  been  reported  in 
macrophages  and  T-lymphocytes  as  well  as  in  mammary  epithelial  cells,  is  dramatically  upregulated  at 
day  12  of  pregnancy  and  remains  high  through  day  7  of  involution  (Fig.  5).  Though  declining  somewhat 
by  day  28  of  postlactational  involution,  Eta-1  expression  remains  markedly  elevated  compared  to  age- 
matched  nulliparous  animals.  In  aggregate,  our  data  suggest  that  parity  induces  persistent  increases  in 
populations  of  hematopoietic  cells  present  within  the  mammary  gland,  as  well  as  changes  in  cytokine 
expression  within  the  mammary  epithelium  itself. 


Task  2:  Months  6-24:  Determine  the  spatial  pattern  of  expression  ofcDNA  markers  as  a  function  of 

reproductive  endocrine  history. 

This  task  has  been  completed  on  schedule.  The  parity-specific  expression  of  a  given  gene  could 
reflect  a  global  increase  in  its  expression  in  all  expressing  cells,  an  increase  in  the  percentage  of 
expressing  cells  in  the  breast,  or  both.  In  order  to  distinguish  between  these  possibilities,  we  are 
performing  in  situ  hybridization  for  those  genes  that  have  been  identified  as  being  differentially 
expressed  as  a  function  of  parity.  As  an  example,  the  spatial  expression  of 
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Fig.  6;  Spatial  expression  of  amphiregulin  during  development 

amphiregulin  is  shown  during  development  (Fig.  6).  This  analysis  demonstrates  that  the  expression  of 
amphiregulin  decreases  in  each  expressing  cell,  and  also  that  the  number  of  amphiregulin-expressing 
cells  decreases  as  a  result  of  parity. 
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Similar  approaches  have  been  applied  to  the  other  differentially  expressed  genes  isolated  in  this 
panel.  For  example,  analysis  of  the  spatio-temporal  pattern  of  gene  expression  for  the  hematopoietic 
markers  described  above  has  also  been  determined  (Fig.  7).  These  finding  clearly  demonstrate  that 
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Fig.  7:  Spatial  expression  of  hematopoietic  gene  expression  during  development 


parity-dependent  increases  in  expression  of  MME,  k-LC  and  Eta-1  are  due  to  increases  in  the  number  of 
cells  that  express  these  genes. 

Task  3:  Months  6-24:  Determine  the  relationship  between  the  expression  of  cDNA  markers  and 

reproductive  endocrine  variables  known  to  affect  breast  cancer  risk. 

This  task  has  been  completed  on  schedule.  We  have  determined  how  specific  parity-related 
reproductive  variables  alter  the  expression  of  selected  genes  in  the  murine  breast  that  we  have  identified 
as  being  differentially  expressed  as  a  function  of  parity.  For  example,  epidemiologic  findings  suggest 
that  those  parity-induced  changes  in  the  mammary  gland  that  are  responsible  for  proteetion  against 
breast  caneer  are  likely  to  be  permanent.  Aecordingly,  ehanges  in  gene  expression  that  are  involved  in 
this  protective  effect  would  be  predicted  to  persist  for  periods  of  involution  greater  than  4  weeks.  To 
determine  whether  the  parity-dependent  changes  in  gene  expression  identified  in  this  study  persist  for 
longer  period  of  involution,  we  analyzed  cohorts  of  mice  that  were  mated  at  4  weeks  of  age,  and  then 
underwent  21  days  of  laetation  and  either  4,  16,  or  30  weeks  of  postlactational  involution.  Northern 
analysis  revealed  that  expression  levels  of  lactoferrin,  k-LC,  TGF-fi3,  clusterin,  and  Eta-1,  were  all 
consistently  upregulated  in  the  mammary  glands  of  parous  animals  compared  to  age-matched 
nulliparous  controls  for  up  to  30  weeks  of  postlactational  involution  (Fig.  8).  These  findings  indicate 
that  for  at  least  a  subset  of  the  genes  identified  in  this  study  parity-induced  ehanges  in  gene  expression 
are  essentially  permanent.  In  addition,  we  have  assessed  the  impact  of  time  of  first  full-term  pregnaney, 
multiparity,  duration  of  lactation,  and  age,  on  the  expression  of  the  genes  isolated  in  these  studies. 
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Fig.  8:  Expression  of  parity-specific  genes  as  a  function  of  duration  of  postlactational  regression 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Identification  of  genes  that  are  down-regulated  as  a  function  of  parity. 

•  Identification  of  genes  that  are  up-regulated  as  a  function  of  parity. 

•  Analysis  of  developmental  and  spatial  patterns  of  gene  expression  for  genes  regulated  by  parity. 

•  Successful  application  of  oligonucleotide  microarray  technology  to  gene  expression  in  whole  organs. 

•  Validation  of  the  overall  approach  proposed  in  this  application. 
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REPORTABLE  OUTCOMES 

•  A  manuscript  describing  the  results  of  this  study  has  been  submitted  for  publication.  This 
manuscript  is  included  in  the  Appendix  to  this  report. 


CONCLUSIONS 

Each  of  the  technical  objectives  and  tasks  proposed  in  this  project  have  been  completed  on  schedule. 
Our  principal  advance  has  been  the  identification  of  molecular  changes  that  occur  in  the  mammary  gland 
as  a  consequence  of  parity.  Such  genetic  changes  have  not  previously  been  identified.  We  have  used 
high-density  oligonucleotide  microarrays  to  analyze  the  impact  of  early  first  full-term  pregnancy  on 
global  gene  expression  profiles  within  the  murine  mammary  gland.  This  approach  has  led  to  the 
identification  of  a  panel  of  genes  whose  expression  in  the  mammary  gland  is  persistently  altered  as  a 
consequence  of  parity  in  multiple  strains  of  mice  as  well  as  in  two  widely  used  rat  models  for  parity- 
induced  protection  against  breast  cancer.  Our  findings  demonstrate  that  parity  induces  the  persistent 
downregulation  of  multiple  genes  encoding  epithelial  growth  factors  as  well  as  the  persistent 
upregulation  of  the  growth-inhibitory  molecule,  TGF-fi3,  and  several  of  its  downstream  targets.  Our 
findings  further  indicate  that  parity  results  in  a  persistent  increase  in  the  differentiated  state  of  the 
mammary  gland  as  well  as  permanent  changes  in  the  hematopoietic  cell  types  resident  within  the  gland. 
These  findings  define  a  developmental  state  of  the  mammary  gland  that  is  refractory  to  carcinogenesis 
and  suggest  novel  hypotheses  for  mechanisms  by  which  parity  may  modulate  breast  cancer  risk.  We 
believe  that  these  studies  represent  an  important  step  towards  determining  the  mechanisms  by  which 
breast  cancer  susceptibility  is  modulated  by  reproductive  history.  This  is  particularly  important  since 
the  identification  and  use  of  intermediate  molecular  endpoints  that  accurately  identify  changes  in  the 
breast  associated  with  changes  in  breast  cancer  risk  will  ultimately  facilitate  the  development  of 
chemopreventative  regimens  that  mimic  parity-induced  protection  against  breast  cancer. 
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ABSTRACT 

Epidemiologic  studies  have  repeatedly  demonstrated  that  women  who  undergo  a  first  full- 
term  pregnancy  early  in  life  have  a  significantly  reduced  lifetime  risk  of  breast  cancer.  Similar 
to  humans,  rodents  that  have  previously  undergone  a  full-term  pregnancy  are  highly  resistant  to 
carcinogen-induced  breast  cancer  compared  to  age-matched  nulliparous  controls.  Unfortunately, 
little  progress  has  been  made  over  the  past  30  years  towards  understanding  the  cellular  and 
molecular  basis  of  parity-induced  protection  against  breast  cancer.  We  have  used  high-density 
DNA  oligonucleotide  microarrays  to  identify  genes  whose  expression  in  the  murine  mammary 
gland  is  persistently  altered  as  a  consequence  of  an  early  first  full-term  pregnancy.  The  panel  of 
differentially  expressed  genes  that  we  have  isolated  in  FVB  mice  reproducibly  distinguishes 
between  the  nulliparous  and  parous  states  of  the  mammary  gland  in  multiple  strains  of  mice  and 
rats.  We  find  that  parity  results  in  the  persistent  down-regulation  of  multiple  genes  encoding 
epithelial  growth  factors  as  well  as  the  persistent  upregulation  of  the  growth-inhibitory  molecule, 
TGF-fi3,  and  several  of  its  transcriptional  targets.  Our  studies  further  indicate  that  parity  results 
in  a  persistent  increase  in  the  differentiated  state  of  the  mammary  gland  as  well  as  persistent 
changes  in  the  immunologic  environment  and  hematopoietic  cell  types  resident  within  the  gland. 
These  data  provide  the  first  molecular  description  of  a  developmental  state  of  the  mammary 
gland  that  is  associated  with  reduced  cancer  risk  and  suggest  novel  hypotheses  for  the 
mechanistic  basis  by  which  parity  may  modulate  breast  cancer  risk. 
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INTRODUCTION 

Numerous  epidemiologic  studies  have  shown  that  women  who  undergo  a  full-term 
pregnancy  early  in  life  have  a  significantly  reduced  lifetime  risk  of  breast  cancer  (1-10).  In 
contrast,  women  who  undergo  their  first  full-term  pregnancy  after  age  30  have  a  risk  of  breast 
cancer  that  is  actually  higher  than  that  of  nulliparous  women.  This  suggests  that  parity-induced 
protection  against  breast  cancer  is  principally  dependent  upon  the  timing  of  a  first  full-term 
pregnancy,  rather  than  on  its  occurrence  per  se.  Other  reproductive  variables  such  as  multiparity 
and  lactation  have  also  been  shown  to  reduce  breast  cancer  risk,  but  in  a  manner  that  is 
independent  of  age  at  first  childbirth.  Moreover,  these  effects  are  generally  modest  in 
comparison  to  age  at  first  full-term  pregnancy  (1,  6,  11-16). 

Interestingly,  women  from  different  countries  and  ethnic  groups  exhibit  a  similar  degree  of 
parity-induced  protection  against  breast  cancer  regardless  of  whether  the  regional  incidence  of 
this  malignancy  is  high,  as  in  Western  countries,  or  low  as  in  the  Far  East.  This  suggests  that  the 
reduction  in  breast  cancer  risk  associated  with  early  first  full-term  pregnancy  does  not  result 
from  extrinsic  factors  specific  to  a  particular  environmental,  genetic,  or  socioeconomic  setting, 
but  rather  from  an  intrinsic  effect  of  parity  on  the  biology  of  the  breast.  In  principle,  this 
protective  effect  could  result  from  the  pregnancy-driven  terminal  differentiation  of  a 
subpopulation  of  target  cells  at  increased  risk  for  carcinogenesis,  from  the  preferential  loss  of 
target  cells  during  postlaetational  involution,  or  from  a  permanent  endocrine  change  that 
indirectly  decreases  breast  cancer  risk  by  altering  either  the  hormonal  environment  or  the 
hormonal  responsiveness  of  cells  in  the  mammary  gland  [Russo,  1978  #1 181;  Russo,  1979  #719; 
Russo,  1982  #1182;  Thordarson,  1995  #1748;  Swanson,  1995  #1739;  (17-19)].  To  date, 
however,  little  evidence  exists  to  support  of  any  of  these  hypotheses  at  the  cellular  or  molecular 
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level.  As  such,  more  than  30  years  after  the  landmark  description  of  this  phenomenon,  the 
mechanisms  that  underlie  parity-induced  protection  against  breast  cancer  remain  largely 
unknown. 

Like  humans,  both  rats  and  mice  exhibit  parity-induced  protection  against  breast  cancer. 

For  example,  administration  of  the  carcinogens  7,12-dimethylbenz(a)anthracene  or  N- 
methylnitrosourea  to  nulliparous  rats  induces  mammary  adenocarcinomas  that  are  hormone- 
dependent  and  histologically  similar  to  human  breast  tumors  [Russo,  1978  #1181;  Russo,  1979 
#719;  Russo,  1977  #339;  Moon,  1969  #261;  Moon,  1981  #262;  Huggins,  1959  #174;  Huggins, 
1959  #170;  (20-22)Guzman,  1999  #1558;  (23)Yang,  1999  #1786].  In  contrast,  rats  that  have 
previously  undergone  a  full-term  pregnancy  are  highly  resistant  to  the  induction  of  breast  cancer 
by  carcinogen  administration.  Similar  to  rats,  mice  that  have  undergone  an  early  first  full-term 
pregnancy  have  also  been  shown  to  be  less  susceptible  to  carcinogen-induced  mammary  tumors 
than  age-matched  nulliparous  controls  (24).  These  studies  indicate  that  key  epidemiological 
features  of  the  influence  of  reproductive  history  on  breast  cancer  risk  in  humans  are  mirrored  in 
rodent  model  systems.  The  use  of  animal  models  to  study  parity-induced  protection  against 
breast  cancer  is  further  facilitated  by  the  many  similarities  in  structure,  function,  and 
development  that  exist  between  the  human  and  rodent  mammary  gland  (25,  26).  Since  the 
mechanisms  underlying  parity-induced  changes  in  breast  cancer  risk  are  likely  to  involve 
complex  genetic  and  epigenetic  events,  animal  models  that  recapitulate  epidemiological  findings, 
permit  critical  aspects  of  reproductive  history  to  be  rigorously  controlled,  reduce  genetic 
variation,  and  permit  the  examination  of  molecular  and  cellular  events  at  defined  developmental 
stages  of  interest  in  normal  tissue  will  be  critical  for  understanding  this  phenomenon. 

Despite  long-standing  evidence  for  the  differential  susceptibility  of  the  parous  and 
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nulliparous  breast  to  carcinogenesis,  a  comprehensive  analysis  of  the  molecular  and  cellular 
changes  induced  in  the  breast  by  parity  has  not  been  previously  reported.  Such  information 
would  not  only  define  a  protected  state  of  the  mammary  gland  at  the  molecular  level,  but  could 
also  provide  insight  into  the  pathways  that  underlie  parity-induced  protection.  In  addition, 
identifying  a  panel  of  molecular  markers  whose  expression  is  reproducibly  altered  by  early  first 
full-term  pregnancy  would  provide  candidate  intermediate  molecular  endpoints  by  which  to 
monitor  the  efficacy  of  pharmacological  interventions  designed  to  mimic  this  naturally  occurring 
protective  event. 

In  this  report  we  have  used  high-density  oligonucleotide  microarrays  to  analyze  the 
impact  of  early  first  full-term  pregnancy  on  global  gene  expression  profiles  within  the  murine 
mammary  gland.  This  approach  has  led  to  the  identification  of  a  panel  of  genes  whose 
expression  in  the  mammary  gland  is  persistently  altered  as  a  consequence  of  parity  in  multiple 
strains  of  mice  as  well  as  in  two  widely  used  rat  models  for  parity-induced  protection  against 
breast  cancer.  Our  findings  demonstrate  that  parity  induces  the  persistent  downregulation  of 
multiple  genes  encoding  epithelial  growth  factors  as  well  as  the  persistent  upregulation  of  the 
growth-inhibitory  molecule,  TGF-fi3,  and  several  of  its  downstream  targets.  Our  findings  further 
indicate  that  parity  results  in  a  persistent  increase  in  the  differentiated  state  of  the  mammary 
gland  as  well  as  permanent  changes  in  the  hematopoietic  cell  types  resident  within  the  gland. 
These  findings  define  a  developmental  state  of  the  mammary  gland  that  is  refractory  to 
carcinogenesis  and  suggest  novel  hypotheses  for  mechanisms  by  which  parity  may  modulate 
breast  cancer  risk. 
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RESULTS 

Parity  Results  in  Permanent  Morphological  Changes  in  the  Mammary  Gland 

In  humans,  parity-induced  changes  in  breast  cancer  susceptibility  are  accompanied  by 
morphological  alterations  in  the  mammary  gland  (26-28).  These  structural  changes  persist 
throughout  life  and  have  been  interpreted  as  a  parity-induced  increase  in  the  differentiated  state 
of  the  breast  (26-28).  Parity-induced  changes  in  morphology  and  cancer  susceptibility  have  also 
been  observed  the  mammary  glands  of  rats  and  mice  [Sinha,  1988  #1903;  Sivaraman,  1998 
#1724;  Russo,  1980  #1703;  Moon,  1969  #261(22,  24)].  To  extend  these  findings,  we  compared 
the  morphology  of  mammary  glands  from  nulliparous  and  parous  Sprague-Dawley  rats,  a  widely 
used  model  for  parity-induced  protection  against  breast  cancer,  with  those  from  C57B1/6  and 
FVB  mice. 

Following  mating  at  an  age  corresponding  to  the  onset  of  puberty,  rats  and  mice  were 
allowed  to  undergo  a  single  round  of  pregnancy,  21  days  of  lactation,  and  28  days  of 
postlactational  involution.  Examination  of  carmine-stained  whole  mounts  prepared  from  these 
animals  demonstrated  that  the  epithelial  trees  of  both  the  rat  and  mouse  parous  involuted 
mammary  gland  are  more  highly  branched  than  those  of  age-matched  nulliparous  littermates 
(Fig.  1).  These  differences  persist  for  at  least  30  weeks  of  postlactational  involution  indicating 
that  the  morphological  changes  induced  by  parity  are  essentially  permanent  (data  not  shown). 
Despite  the  marked  differences  in  ductal  branching  patterns  between  nulliparous  and  parous 
mammary  glands,  hematoxylin  and  eosin-stained  sections  reveal  that  the  relative  amount  of 
epithelium  contained  within  the  mammary  fat  pad  is  similar  in  both  developmental  states  (Fig.  1 
and  data  not  shown).  Thus,  architectural  differences  between  the  nulliparous  and  parous  mouse 
mammary  gland  are  easily  distinguishable,  occur  in  multiple  strains,  are  conserved  among  rodent 
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species,  and  are  analogous  to  those  that  have  been  described  in  the  human  breast  (26,  27). 

Microarray  Analysis  of  Parity-Induced  Changes  in  Gene  Expression 

Since  parity-induced  changes  in  both  the  structure  of  the  mammary  gland  and  its 
susceptibility  to  cancer  are  conserved  between  humans  and  rodents,  we  sought  to  define  the 
molecular  differences  between  the  nulliparous  and  parous  mammary  gland  using  high  density 
oligonucleotide  microarrays.  These  studies  were  intended  to  define  a  developmentally  protected 
state  of  the  mammary  gland  at  the  molecular  level  as  well  as  facilitate  the  identification  of  parity- 
induced  changes  in  the  abundance  of  different  cell  types  within  the  mammary  gland.  While  we 
anticipated  that  many  of  the  differentially  expressed  genes  identified  by  this  approach  might  not 
play  a  causal  role  in  parity-induced  protection,  we  nevertheless  considered  it  likely  that 
differentially  expressed  genes  would  provide  insight  into  parity-induced  alterations  in  the 
mammary  gland,  including  those  that  are  responsible  for  the  resistance  of  the  parous  gland  to 
carcinogenesis. 

Oligonucleotide  microarray  expression  profiling  was  performed  in  triplicate  on  pooled 
mammary  gland  samples,  each  of  which  was  derived  from  15-20  animals  to  control  for  sources 
of  biological  variation  (29).  Since  stromal-epithelial  interactions  have  been  clearly  shown  to 
affect  the  behavior  of  mammary  epithelial  cells,  expression  changes  were  profiled  in  intact 
mammary  glands  with  the  exception  that  the  lymph  node  present  in  the  number  4  mammary 
gland  was  removed  [Wiesen,  1999  #1775(30-33)].  Although  the  morphological  changes 
characteristic  of  postlactational  involution  are  essentially  complete  after  14  days,  we  chose  to 
profile  parous  animals  after  28  days  of  involution  to  facilitate  the  identification  of  persistent 
changes  in  gene  expression  due  to  parity  rather  than  acute  changes  in  gene  expression  due  to  the 
processes  of  pregnancy,  lactation  or  involution  per  se. 
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Three  independent,  age-matched  pools  of  total  RNA  derived  from  nulliparous  (15  wk 
GOTO)  and  parous  (15  wk  GlPl)  cohorts  were  hybridized  to  high-density  oligonucleotide 
microarrays  representing  approximately  5,300  murine  genes  and  ESTs  (Mu6500).  Affymetrix 
comparative  algorithms  were  used  to  identify  genes  whose  expression  levels  consistently 
changed  as  a  consequence  of  parity.  Northern  hybridization  analysis  performed  on  independent 
nulliparous  and  parous  mammary  gland  samples  confirmed  the  differential  expression  of  14/14 
genes  identified  by  Affymetrix  algorithms  as  being  differentially  expressed  in  each  of  the  three 
independent  microarray  comparisons  (Fig.  2  and  data  not  shown).  This  approach,  combined  with 
Northern  analysis  of  candidate  genes  identified  by  microarray  analysis  as  being  differentially 
expressed  in  two  out  of  three  sets  of  nulliparous  and  parous  samples,  resulted  in  the  identification 
of  41  differentially  expressed  genes  (Table  1  and  Fig.  2).  These  include  10  genes  that  are 
preferentially  expressed  in  the  nulliparous  mammary  gland  as  well  as  3 1  genes  that  are 
preferentially  expressed  in  the  parous  mammary  gland.  Notably,  expression  levels  of  the 
epithelial  cell  marker,  cytokemtinlS,  do  not  differ  between  the  nulliparous  and  parous  mammary 
gland,  indicating  that  the  observed  changes  in  gene  expression  are  not  merely  a  consequence  of 
the  expansion  or  contraction  of  the  epithelial  cell  compartment  (Fig.  2). 

Differentially  Expressed  Genes  Distinguish  Parous  and  Nulliparous  Tissues  in  Mice  and 
Rats 

Having  identified  a  panel  of  genes  that  were  differentially  expressed  in  a  parity- 
dependent  manner  in  an  index  group  of  FVB  mice,  we  asked  whether  the  expression  patterns  of 
these  genes  were  sufficient  to  consistently  distinguish  nulliparous  and  parous  mammary  tissues 
harvested  from  independent  groups  of  FVB  mice  as  well  as  from  additional  strains  of  mice.  To 
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this  end,  we  generated  12  additional  independent  sets  of  pooled  mammary  gland  RNA  samples 
from  parous  and  nulliparous  FVB  and  Balb/c  mice.  These  samples  were  subjected  to  high- 
density  oligonucleotide  microarray  analysis  on  MGU74A  mouse  arrays  and  clustered  based  on 
expression  profiles  for  the  differentially  expressed  genes  identified  in  our  original  Mu6500  array 
analysis.  Cluster  analysis  performed  in  a  blinded  manner  demonstrated  that  the  expression 
patterns  for  the  41  genes  identified  in  this  study  were  sufficient  to  correctly  distinguish 
nulliparous  and  parous  mammary  gland  samples  harvested  from  independent  sets  of  FVB  mice 
(Fig.  3A).  Furthermore,  the  expression  patterns  for  these  same  genes  were  also  sufficient  to 
correctly  distinguish  nulliparous  and  parous  mammary  gland  samples  harvested  from  Balb/c 
mice  (Fig.  3B).  Given  that  different  generations  of  microarrays  composed  of  different  probe  sets 
were  used  to  perform  these  analyses,  our  findings  suggest  that  this  panel  of  genes  accurately  and 
reproducibly  identifies  parity-induced  changes  in  the  mouse  mammary  gland  (Fig.  3A). 

Since  Sprague-Dawley  and  Lewis  rats  represent  the  most  widely  used  models  of  parity- 
induced  protection  against  breast  cancer,  we  were  further  interested  in  determining  whether 
genes  identified  as  having  a  parity-dependent  pattern  of  expression  in  FVB  mice  would  be 
sufficient  to  predict  correctly  the  reproductive  histories  of  rats.  As  such,  6  independent  sets  of 
pooled  mammary  gland  RNA  samples  from  parous  and  nulliparous  Lewis  rats  were  analyzed  on 
RGU34A  high-density  rat  microarrays  To  facilitate  comparison  of  mouse  data  sets  genes  to 
microarray  expression  data  obtained  from  rats,  genes  identified  as  being  expressed  in  parity- 
dependent  maimer  in  FVB  mice  were  mapped  via  Homologene  to  the  rat  genome.  Cluster 
analysis  of  rat  microarray  data  revealed  that  the  panel  of  genes  identified  on  the  basis  of  their 
parity-dependent  expression  in  FVB  mice  was  sufficient  to  correctly  distinguish  mammary  gland 
samples  harvested  from  nulliparous  and  parous  Lewis  rats  (Fig.  3C).  These  findings  demonstrate 
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that  the  expression  patterns  of  the  genes  that  we  have  isolated  are  reproducibly  and  persistently 
altered  as  a  consequence  of  parity,  and  that  these  parity-induced  changes  in  gene  expression  are 
conserved  in  different  mouse  strains  and  rodent  species  that  exhibit  parity-induced  protection 
against  breast  cancer. 

Functional  Gene  Categories  Accurately  Prediet  Reproductive  History 

Examination  of  the  genes  identified  as  being  differentially  expressed  as  a  consequence  of 
parity  revealed  several  distinct  functional  categories  (Table  1).  These  include:  growth-promoting 
molecules  such  as  amphiregulin,  pleiotrophin,  and  insulin-like  growth  factor  I;  molecules  related 
to  epithelial  differentiation  such  as  milk  proteins;  molecules  expressed  by  hematopoietic  cells 
such  as  B-cells,  T-cells  and  macrophages;  and  molecules  involved  in  the  TGF-B  pathway.  In 
light  of  this  observation,  we  asked  whether  smaller  subsets  of  genes  representing  these  four 
functional  gene  categories  (growth  factor-related:  4  genes;  differentiation-related:  7  genes; 
immune-related:  15  genes;  and  TGF-B-related:  4  genes)  would  be  sufficient  to  correctly 
distinguish  nulliparous  and  parous  mammary  gland  samples  in  a  blinded  manner. 

Strikingly,  even  when  considered  in  isolation,  gene  expression  patterns  within  each  of 
these  four  subgroups  were  found  to  be  sufficient  to  identify  correctly  parous  and  nulliparous 
mammary  samples  derived  from  FVB  mice,  Balb/c  mice  and  Fewis  rats  (Fig.  3D,  3E  and  data 
not  shown).  These  findings  demonstrate  that  expression  changes  within  each  of  these  four 
functional  gene  categories  robustly  and  independently  distinguish  between  the  parous  and  the 
nulliparous  states  of  the  mammary  gland  in  different  rodent  strains  and  species  that  exhibit 
parity-induced  protection  against  breast  cancer.  Our  observations  further  suggest  that  the 
downregulation  of  specific  genes  involved  in  epithelial  proliferation,  and  the  upregulation  of 
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genes  involved  in  epithelial  differentiation,  immune  regulation,  and  TGF-B-mediated  grovvth 
inhibition  represent  cardinal  features  of  parity-induced  changes  in  the  mammary  gland. 

Differentiation  Markers  are  Preferentially  Expressed  in  the  Parous  Mammary  Gland 

A  prominent  category  of  genes  whose  expression  was  elevated  in  the  parous  involuted 
mammary  gland  included  markers  for  mammary  epithelial  differentiation  such  as  a-casein,  B- 
casein,  y-casein,  K-casein,  whey  acidic  protein,  lactoferrin  and  a-lactalbumin  (Table  1).  Cluster 
analysis  of  microarray  data  obtained  for  6  pooled  mammary  samples  harvested  from  independent 
sets  of  FVB  mice  demonstrated  that  the  expression  patterns  of  these  seven  genes  were  sufficient 
to  correctly  distinguish  nulliparous  and  parous  FVB  samples  in  a  blinded  manner  (Fig.  3  and 
Table  1).  The  differential  expression  patterns  of  these  7  genes  were  also  sufficient  to  distinguish 
correctly  nulliparous  and  parous  mammary  samples  harvested  from  Balb/c  mice  as  well  as  from 
Lewis  rats  (Fig.  3  and  data  not  shown).  Additional  genes  associated  with  epithelial 
differentiation,  such  as  connexin  26,  were  also  found  to  be  persistently  upregulated  by  parity 
(Table  1).  These  observations  indicate  that  the  parity-dependent  upregulation  of  markers  for 
mammary  epithelial  differentiation  is  conserved  among  different  mouse  strains  and  rodent 
species  that  exhibit  parity-induced  protection  against  breast  cancer.  These  findings  provide  the 
first  molecular  evidence  supporting  the  hypothesis  that  parity  results  in  a  persistent  increase  in 
the  differentiated  state  of  the  mammary  epithelium  (34).  In  addition,  the  preferential  expression 
of  ADFP,  which  is  upregulated  in  differentiated  adipocytes,  in  the  parous  involuted  gland 
suggests  that  stromal  compartments  of  the  mammary  gland  may  also  become  more  differentiated 
as  a  consequence  of  parity  (35). 
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Parity-Induced  Changes  in  Hematopoietic  Cells  in  the  Mammary  Gland 

A  second  category  of  genes  whose  expression  was  elevated  in  the  parous  involuted 
mammary  gland  included  those  that  are  specifically  expressed  in  hematopoietic  cells  such  as  B- 
lymphocytes  {kappa  light  chain  and  the  IgG,  IgM,  and  IgA  heavy  chains),  T-lymphocytes  {T-cell 
death  associated  gene),  and  macrophages  {macrophage  metalloelastase  and  macrophage 
expressed  gene  1)  (Table  1  and  Fig.  4A).  Additional  genes,  including  Early  T-cell  activation 
protein  {Eta-1),  LPS-Binding  protein  and  Lipocalin-2,  that  are  either  expressed  by  hematopoietic 
cells  or  are  chemoattractants  for  these  cell  types  were  also  persistently  upregulated  as  a 
consequence  of  parity  (36-39).  Similar  to  markers  of  epithelial  differentiation,  the  expression 
patterns  of  this  subset  of  immune-related  genes  were  able  in  a  blinded  manner  to  identify 
correctly  parous  and  nulliparous  mammary  gland  samples  in  FVB  mice,  Balb/c  mice  and  Lewis 
rats  (Fig.  3  and  data  not  shown).  The  marked  upregulation  of  these  genes  in  the  parous  gland 
suggests  that  cells  of  B-lymphocyte,  T-lymphocyte,  and  macrophage  lineages  may  be  more 
abundant  in  the  parous  mammary  gland. 

To  investigate  this  hypothesis  further,  we  examined  the  expression  patterns  for  several 
immune-related  genes  during  stages  of  postnatal  mammary  development  representing  puberty, 
pregnancy,  lactation  and  involution.  Northern  analysis  revealed  that  expression  of  the  B-cell- 
specific  gene,  K-light  chain  {k-LC),  is  first  detected  during  lactation  with  elevated  levels  of 
expression  persisting  throughout  postlactational  involution  (Fig.  4B).  The  increase  in  k-LC 
expression  during  lactation  is  partially  masked  by  the  dilutional  effects  that  result  from  the 
massive  increase  in  milk  protein  gene  expression  that  occurs  during  this  period,  as  evidenced  by 
the  apparent  decrease  in  _-actin  expression  that  occurs  during  this  same  period  (Fig.  4B)  (40, 

41).  In  situ  hybridization  confirmed  the  upregulation  of  k-LC  expression  during  lactation  and 
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further  revealed  that  this  upregulation  is  due  to  an  increase  in  the  number  of  k-LC  expressing 
cells  (Fig.  4C  and  data  not  shown).  The  spatio-temporal  pattern  of  k-LC  expression  is  consistent 
with  the  reported  influx  of  lymphocytes  into  the  breast  that  occurs  during  lactation  (42-44). 
Surprisingly,  however,  this  lymphocyte  population  appears  to  persist  in  the  fiilly  involuted  gland. 

Macrophage  metalloelastase  (MME  or  MMP-12)  was  also  found  to  be  persistently 
upregulated  in  the  mammary  gland  as  a  consequence  of  parity  (Fig.  4A).  MME  is  a  secreted 
metalloprotease  that  cleaves  plasminogen  to  generate  angiostatin,  a  potent  inhibitor  of 
endothelial  cell  proliferation  (45-47).  Northern  hybridization  demonstrated  that  MME 
expression  levels  increase  dramatically  at  day  7  of  postlactational  involution  and  remain  elevated 
compared  to  age-matched  nulliparous  controls  following  28  days  of  involution  (Fig.  4B).  In  situ 
hybridization  revealed  that  MME  expression  is  restricted  to  isolated  cells  within  the  mammary 
stroma  at  day  7  of  involution  (Fig.  4C  and  data  not  shown).  This  finding  is  consistent  with  this 
gene's  reported  expression  in  macrophages  and  with  previous  evidence  that  macrophages  are 
recruited  to  the  breast  during  involution  where  they  participate  in  the  clearance  of  post-apoptotic 
debris  (48,  49).  Interestingly,  by  day  28  of  involution  foci  of  MME  expression  became  tightly 
associated  with  the  epithelial  compartment  reflecting  either  a  persistent  population  of 
macrophages  residing  within  the  parous  epithelium  or  a  subset  of  epithelial  cells  expressing 
MME  (Fig.  AC). 

Finally,  we  examined  the  developmental  basis  for  the  persistent  parity-dependent 
upregulation  of  Eta-1  expression  (Fig.  4A).  Northern  analysis  revealed  that  Eta-1  expression, 
which  has  been  reported  in  macrophages  and  T-lymphocytes  as  well  as  in  mammary  epithelial 
cells,  is  dramatically  upregulated  at  day  12  of  pregnancy  and  remains  high  through  day  7  of 
involution  (Fig.  4B)  (50,  51).  Though  declining  somewhat  by  day  28  of  postlactational 
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involution,  Eta-1  expression  remains  markedly  elevated  compared  to  age-matched  nulliparous 
animals.  Notably,  in  situ  hybridization  analysis  demonstrated  that  Eta-1  expression  in  the  parous 
involuted  gland  is  restricted  to  a  subset  of  mammary  epithelial  cells  (Fig.  4C  and  data  not 
shown).  In  aggregate,  our  data  suggest  that  parity  induces  persistent  increases  in  populations  of 
hematopoietic  cells  present  within  the  mammary  gland,  as  well  as  changes  in  cytokine 
expression  within  the  mammary  epithelium  itself. 

Parity  Results  in  a  Decrease  in  Growth  Factor  Expression 

Interestingly,  genes  that  encode  growth  regulatory  molecules  constitutex  more  than  half 
of  genes  that  were  found  to  be  persistently  downregulated  by  parity  (Table  1).  These  include 
amphiregulin  (Areg),  pleiotrophin  (Ptn),  insulin-like  growth  factor  1  (Igfl),  leptin  (Ob),  and 
thyroid  stimulating  hormone  receptor  (TshR)  (Table  1  and  Fig.  5  A).  Cluster  analysis  further 
demonstrated  that  the  expression  patterns  of  these  5  genes  were  sufficient  in  a  blinded  manner  to 
identify  correctly  parous  and  nulliparous  mammary  samples  in  FVB  mice,  Balb/c  mice  and 
Lewis  rats  (Fig.  3).  This  indicates  that  downregulation  of  the  genes  encoding  these  growth 
factors  is  a  characteristic  feature  of  parity-induced  changes  in  the  rodent  mammary  gland. 

We  investigated  the  developmental  expression  pattern  of  the  epidermal  growth  factor 
receptor  ligand,  Areg,  and  the  heparin-binding  mitogen,  Ptn,  to  determine  the  basis  for  their 
preferential  expression  in  the  nulliparous  gland.  Strikingly,  the  temporal  patterns  of  expression 
for  these  molecules  are  virtually  identical  during  postnatal  mammary  development  (Fig.  5B). 
Steady  state  levels  of  Areg  and  Ptn  mRNA  are  dramatically  upregulated  in  the  female  mammary 
gland  between  2  weeks  and  5  weeks  of  age,  a  period  corresponding  to  the  onset  of  ductal 
morphogenesis  (David  Lee  ref).  Levels  of  Areg  and  Ptn  expression  remain  relatively  constant 
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throughout  nulliparous  development  and  early  pregnancy,  then  decrease  sharply  by  mid¬ 
pregnancy  (day  12)  and  remain  low  throughout  lactation  and  involution  (Fig.  5B).  In  situ 
hybridization  analysis  confirmed  the  pregnancy-induced  downregulation  of  Areg  expression  as 
well  as  the  persistent  downregulation  of  Areg  expression  throughout  the  epithelial  compartment 
of  the  parous  gland  (Fig.  5C).  Together  with  the  observed  decreases  in  Ptn  and  Igfl  expression 
levels,  our  data  suggest  the  possibility  that  parity  results  in  the  downregulation  of  multiple 
pathways  that  stimulate  epithelial  proliferation. 

Parity  Results  in  Increased  Mammary  Expression  of  TGF-B3  and  its  Targets 

Our  microarray  data  indicated  that  steady-state  mRNA  levels  for  the  growth-inhibitory 
cytokine,  TGF-fi3,  were  persistently  upregulated  in  the  mammary  glands  of  FVB  mice  as  a 
consequence  of  parity  (Table  1,  Fig.  4,  and  Fig.  6).  In  addition,  clusterin,  Eta-1  and  Id-2,  each  of 
which  has  previously  been  implicated  as  a  downstream  transcriptional  target  of  TGF-fi3,  were 
also  found  to  be  persistently  upregulated  by  parity  (49, 52-56).  Cluster  analysis  demonstrated 
that  the  expression  patterns  of  these  four  genes  are  sufficient  to  correctly  distinguish  parous  and 
nulliparous  mammary  samples  in  FVB  mice,  Balb/c  mice  and  Lewis  rats.  Northern  and  in  situ 
hybridization  analysis  confirmed  the  parity-dependent  upregulation  of  TGF-fi3,  clusterin,  and 
Eta-1  expression  in  the  parous  mammary  gland,  and  further  demonstrated  that  these  genes  are 
primarily  expressed  in  the  epithelial  compartment  (Fig.  6A  and  6C;  Fig.  4A  and  4C).  Northern 
analysis  of  the  developmental  expression  patterns  of  TGF-J33  and  clusterin  demonstrated  that 
both  genes  exhibit  maximal  expression  at  day  2  of  involution  with  a  secondary  peak  of 
expression  at  day  12  of  pregnancy  (Fig.  6B  and  6C).  These  developmental  expression  profiles 
are  consistent  with  a  role  for  these  molecules  in  cell  death  during  mammary  gland  involution  (49, 
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53,  57).  The  coordinate  elevation  in  expression  of  TGF-fi3  along  with  several  of  its 
transcriptional  targets  suggests  that  the  upregulation  of  TGF-J33  mRNA  in  the  parous  mammary 
gland  may  be  accompanied  by  a  bona  fide  increase  in  TGF-fi3  activity. 

Parity-Induced  Changes  in  Gene  Expression  Persist  in  the  Mammary  Gland 

Epidemiologic  findings  suggest  that  those  parity-induced  changes  in  the  mammary  gland 
that  are  responsible  for  protection  against  breast  cancer  are  likely  to  be  permanent.  Accordingly, 
changes  in  gene  expression  that  are  involved  in  this  protective  effect  would  be  predicted  to 
persist  for  periods  of  involution  greater  than  4  weeks.  To  determine  whether  the  parity- 
dependent  changes  in  gene  expression  identified  in  this  study  persist  for  longer  period  of 
involution,  we  analyzed  cohorts  of  mice  that  were  mated  at  4  weeks  of  age,  and  then  underwent 
21  days  of  lactation  and  either  4, 16,  or  30  weeks  of  postlactational  involution.  Northern 
analysis  revealed  that  expression  levels  of  lactoferrin,  _-LC,  TGF-J33,  clusterin,  and  Eta-1,  were 
all  consistently  upregulated  in  the  mammary  glands  of  parous  animals  compared  to  age-matched 
nulliparous  controls  for  up  to  30  weeks  of  postlactational  involution  (Fig.  7A).  These  findings 
indicate  that  for  at  least  a  subset  of  the  genes  identified  in  this  study  parity-induced  changes  in 
gene  expression  are  essentially  permanent. 

Parity-Induced  Changes  in  Gene  Expression  are  Conserved  Across  Species 

Since  the  phenomenon  of  parity-induced  protection  against  breast  cancer  is  conserved 
among  humans,  rats,  and  mice,  we  predicted  that  the  molecular  changes  that  underlie  this  effect 
would  be  conserved.  As  described  above,  microarray  analysis  of  gene  expression  profiles  in  the 
mammary  glands  of  parous  and  nulliparous  FVB  mice,  Balb/c  mice,  and  Lewis  rats  indicated 
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that  parity-dependent  ehanges  in  gene  expression  are  largely  conserved  in  different  strains  of 
mice  and  rats.  To  extend  these  findings,  we  performed  Northern  analysis  on  pools  of  mammary 
gland  total  RNA  isolated  from  nulliparous  and  parous  129SvEv  and  Balb/c  mice,  as  well  as  from 
nulliparous  and  parous  Sprague-Dawley  rats,  another  widely  used  model  for  parity-induced 
protection  against  breast  cancer.  We  found  that  all  ten  genes  examined  in  129SvEv  and  Balb/c 
mouse  strains  were  differentially  expressed  in  a  manner  similar  to  that  observed  in  FVB  mice 
(Fig.  7B  and  data  not  shown).  Likewise,  each  of  the  nine  genes  examined  in  the  rat  exhibited  a 
parity-dependent  differential  pattern  of  expression  identical  to  that  observed  in  the  mouse  (Fig. 
7C  and  data  not  shown).  When  taken  together  with  the  above  microarray  data,  these  findings 
indicate  that  many  of  the  parity-related  molecular  changes  identified  in  this  study  are  conserved 
among  species  that  display  parity-induced  protection  against  mammary  tumorigenesis. 
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DISCUSSION 

We  have  used  DNA  oligonucleotide  microarrays  to  analyze  the  expression  of  -5,300 
genes  and  ESTs  in  order  to  identify  persistent  changes  in  gene  expression  in  the  murine 
mammary  gland  that  are  induced  by  an  early  first  full-term  pregnancy.  Using  this  approach,  we 
have  isolated  a  panel  of  genes  whose  expression  is  persistently  altered  by  a  reproductive  event 
known  to  reduce  breast  cancer  risk.  The  expression  patterns  of  the  genes  isolated  reproducibly 
distinguish  between  the  nulliparous  and  parous  mammary  gland  in  both  rats  and  mice,  as  well  as 
identify  changes  in  the  abundance  of  specific  cell  types  that  occur  in  the  mammary  gland  as  a 
result  of  parity.  Our  findings  demonstrate  at  the  molecular  level  that  parity  results  in  a  persistent 
increase  in  the  differentiated  state  of  the  mammary  epithelium,  as  has  been  previously  suggested 
based  on  morphological  criteria.  In  addition,  our  data  suggest  several  new  hypotheses  for  the 
mechanistic  basis  of  parity-induced  protection  against  breast  cancer,  including  that  parity  may 
decrease  the  susceptibility  of  the  mammary  epithelium  to  malignant  transformation  by 
downregulating  multiple  growth  promoting  pathways,  upregulating  growth-inhibitory  pathways, 
and/or  changing  the  immune  environment  of  the  breast. 

For  the  majority  of  genes  identified  in  this  study,  parity-induced  changes  in  gene 
expression  were  shown  to  be  independent  of  the  length  of  postlactational  involution, 
demonstrating  that  the  altered  patterns  of  expression  that  we  have  identified  are  persistent,  if  not 
permanent.  Furthermore,  we  have  confirmed  that  parity-induced  changes  in  gene  expression  are 
largely  conserved  in  multiple  stains  of  mice  as  well  as  in  two  well-characterized  rat  models  for 
parity-induced  protection  against  breast  cancer.  Beyond  indicating  that  many  of  the  molecular 
changes  that  we  have  defined  are  conserved  among  rodents,  our  findings  suggest  that  these 
changes  may  also  be  relevant  to  parity-induced  changes  in  the  human  breast. 
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One  of  the  most  striking  findings  of  our  microarray  analysis  was  the  parity-dependent 
down-regulation  of  multiple  genes  involved  in  the  regulation  of  cell  growth  and  proliferation. 
Specifically,  the  reduced  expression  of  Areg,  Ptn,  Igfl,  TshR,  and  Ob  suggest  that  multiple 
mitogenic  pathways  may  be  downregulated  in  the  mammary  gland  as  a  consequence  of  parity. 
Notably,  elevated  expression  of  Areg,  Ptn  and  Igfl  have  each  been  implicated  in  the 
pathogenesis  of  human  breast  cancer  (58).  AREG,  a  ligand  for  the  EGER,  has  been  shown  to  be 
overexpressed  in  35-50%  of  primary  human  breast  cancers  (59-61).  Consistent  with  this,  Areg  is 
overexpressed  in  hyperplastic  stages  of  mammary  tumor  development  in  MMTV-PyMT  and 
MT-TGFa  transgenic  mice  (62),  and  is  a  potent  stimulator  of  anchorage-dependent  growth  in 
nontransformed  MCE- 10  cells  (63).  Moreover,  targeted  inactivation  of  Areg  in  mice  causes  a 
marked  delay  in  ductal  elongation  during  puberty  and  inhibition  of  Areg  expression  in 
transformed  mammary  epithelial  cell  lines  by  antisense  approaches  results  in  growth  inhibition 
and  reduced  tumorigenicity  (64-66).  These  findings  indicate  that  Areg  plays  an  important  role  in 
promoting  mammary  epithelial  cell  proliferation  and  predict  that  the  parity-dependent 
downregulation  of  Areg  could  contribute  to  a  decrease  in  the  susceptibility  of  the  parous 
mammary  gland  to  cancer. 

The  heparin-binding  growth  factor,  pleiotrophin,  was  also  found  to  be  down-regulated  in 
the  parous  mammary  gland.  Pleiotrophin  has  been  implicated  in  angiogenesis  and  mammary 
tumor  progression  and  a  majority  of  primary  human  breast  tumors  display  high  levels  of  PTN 
mRNA  expression,  as  do  carcinogen-induced  rat  mammary  tumors  (67,  68).  Moreover,  Ptn 
overexpression  in  NIH  3T3  cells  induces  transformation  and  confers  the  ability  to  form  tumors  in 
nude  mice  (69).  Conversely,  overexpression  of  a  dominant  negative  PTN  mutant  reverses  the 
transformed  phenotype  of  the  human  breast  cancer  cell  line  MDA-MB-231  (70).  PTN  has  also 
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been  reported  to  stimulate  endothelial  cell  growth,  migration,  and  tube  formation.  Furthermorer, 
tumors  formed  from  PJA^-transfected  MCF-7  breast  carcinoma  cells  injected  into  mice  exhibit 
enhanced  growth,  increased  endothelial  proliferation,  and  increased  vascular  density  compared  to 
control  cells  (71-73).  As  such,  like  Areg  the  downregulation  of  Ptn  represents  a  biologically 
plausible  mechanism  that  could  contribute  to  parity-induced  protection  against  breast  cancer. 

It  has  previously  been  reported  that  parity  results  in  the  downregulation  of  the  EGF 
receptor  and  growth  hormone  receptor  in  rats  (22).  Our  finding  that  the  EFGR  ligand,  Areg,  and 
an  EST  that  is  homologous  to  the  growth  hormone  receptor  are  downregulated  by  parity  in  the 
mouse  provide  additional  support  for  this  observation.  However,  our  finding  that  Igfl,  Ptn,  Ob 
and  TshR,  are  also  downregulated  by  parity  raises  the  possibility  that  other  growth  promoting 
pathways  whose  expression  has  not  previously  been  correlated  with  reproduetive  history  may 
also  contribute  to  parity-induced  protection  against  breast  cancer.  This  hypothesis  is  particularly 
intriguing  in  light  of  recent  findings  that  a  strong  positive  correlation  exists  between  circulating 
IGFl  concentrations  and  breast  cancer  risk  among  premenopausal  women  (74).  As  such,  our 
studies  suggest  that  parity  may  result  in  changes  in  both  systemic  and  local  growth  signaling 
pathways  in  the  mammary  gland  that  may  ultimately  contribute  to  the  establishment  of  a  less 
mitogenic  environment. 

In  addition  to  decreases  in  growth  factor  gene  expression,  we  have  also  identified  a 
parity-dependent  increase  in  the  expression  of  TGF-J33  and  several  of  its  downstream 
transcriptional  targets  in  the  mammary  gland.  The  well-described  role  of  TGF-J33  in  growth 
inhibition  makes  the  upregulation  of  this  pathway  a  biologically  plausible  contributing  factor  to 
parity-induced  protection  against  breast  cancer.  While  our  data  do  not  directly  demonstrate 
activation  of  the  TGF-J33  signaling  pathway  in  the  parous  mammary  gland,  the  coordinate 
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regulation  of  downstream  transcriptional  targets  of  TGF-J33  lend  support  this  model.  Further 
investigation  of  the  activation  of  additional  downstream  molecules  in  this  pathway,  such  as 
SMAD  2,  SMAD  3  and  SMAD  4  would  provide  additional  evidence  for  the  upregulation  of  this 
pathway. 

It  has  previously  been  proposed  that  parity-induced  protection  against  breast  cancer  may 
be  mediated  by  an  increased  state  of  differentiation  of  the  parous  mammary  gland  (Russo).  Our 
data  provide  the  first  molecular  evidence  to  support  the  contention  that  the  epithelial 
compartment  of  the  parous  mammary  gland  is  more  differentiated  than  that  of  the  nulliparous 
gland,  as  well  as  evidence  to  suggest  that  parity  may  also  increase  the  differentiated  state  of  the 
stromal  compartment  of  the  mammary  gland.  However,  the  recent  finding  in  rats  that  the 
dopamine  antagonist,  perphenazine,  induces  epithelial  differentiation  yet  does  not  protect  against 
carcinogen-induced  mammary  tumorigenesis  casts  doubt  on  a  central  role  for  differentiation  in 
parity-induced  protection  against  breast  cancer  (18).  As  such,  while  a  persistent  parity-induced 
persistent  increase  in  the  differentiated  state  of  the  mammary  gland  remains  a  biologically 
plausible  mechanism  for  reducing  cancer  susceptibility,  this  finding  may  nevertheless  be 
unrelated  to  the  mechanism  by  which  parity  reduces  breast  cancer  risk.  However,  since  the 
extent  of  epithelial  differentiation  induced  by  perphenazine  may  differ  from  that  induced  by 
parity,  further  studies  will  be  required  to  conclusively  rule  out  differentiation  as  a  contributing 
factor. 

The  finding  that  the  secreted  anti-angiogenesis  factor,  MME,  is  highly  expressed  by  a 
subset  of  epithelial  cells  in  the  involuted  parous  mammary  gland  further  supports  the  hypothesis 
that  parity  induces  changes  in  the  mammary  environment  that  make  it  less  hospitable  for  tumor 
establishment  and  progression  than  the  nulliparous  gland.  The  cleavage  of  plasminogen  by 
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MME  produces  angiostatin,  which  potently  inhibits  endothelial  cell  proliferation  in  vivo  (45-47). 
Gene  transfer  of  MME  into  murine  melanoma  cells  inhibits  angiogenesis  as  well  as  primary 
tumor  growth  (75).  Consistent  with  this  mechanism  of  action,  MME  expression  in  human 
hepatocellular  carcinomas  correlates  highly  with  both  angiostatin  protein  levels  and  patient 
survival  (76).  Thus,  upregulation  of  MME  expression  in  the  parous  mammary  gland  may 
suppress  the  neovascularization  that  is  required  for  tumor  growth. 

Related  to  changes  in  MME  expression,  in  the  course  of  the  above  experiments  we  were 
surprised  to  find  the  parity-induced  upregulation  of  genes  whose  expression  marks  distinct 
classes  of  hematopoietic  cells.  These  observations  suggest  that  parity  results  in  a  persistent 
increase  in  the  number  of  macrophages,  B-lympohcytes  and  T-lymphocytes  that  reside  within  the 
mammary  gland.  In  addition,  we  have  shown  that  parity-induces  changes  in  the  expression  of 
specific  cytokines  in  the  mammary  gland,  at  least  some  of  which  occur  within  the  epithelial 
compartment.  In  particular,  the  dramatic  increase  in  epithelial  expression  of  the  cytokine,  Eta-1, 
which  has  been  shown  to  enhance  migration  of  activated  macrophages,  to  influence  macrophage 
cytokine  production,  and  to  increase  IgG  and  IgM  expression  in  activated  B  cells,  suggests  a 
model  for  paracrine  signaling  mechanisms  by  which  parity-induced  changes  in  epithelial  gene 
expression  may  influence  the  function  of  immune  cells  with  the  gland  (77-79)  (51).  Such 
paracrine  signaling  mechanisms  could  represent  a  potential  route  by  which  parity  might  create  an 
environment  that  is  refractory  to  tumorigenesis.  More  broadly,  however,  our  findings  suggest 
that  parity  may  induce  widespread  changes  in  the  immunological  environment  of  the  mammary 
gland,  potentially  including  tumor  surveillance  mechanisms,  and  may  contribute  to  parity- 
induced  protection  against  breast  cancer. 

Finally,  the  realization  that  specific  reproductive  endocrine  events  alter  breast  cancer  risk 
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in  a  predictable  fashion  raises  the  possibility  that  naturally  occurring  events  known  to  decrease 
breast  cancer  risk  might  be  mimicked  pharmacologically.  The  desire  to  pursue  this  objective  is 
heightened  by  the  fact  that  while  it  is  now  possible  to  identify  women  who  are  at  elevated  risk 
for  developing  breast  cancer,  few  interventions  currently  exist.  As  such,  reducing  breast  cancer 
risk  via  hormonal  manipulations  designed  to  mimic  naturally  occurring  endocrine  events  could 
represent  an  attractive  alternative.  It  is  to  this  end  that  an  early  first  fiill-term  pregnancy  has  been 
proposed  as  a  logical  paradigm  on  which  to  model  the  hormonal  chemoprevention  of  breast 
cancer.  The  achievement  of  this  goal,  however,  has  been  hampered  by  our  lack  of  understanding 
of  the  mechanisms  by  which  reproductive  events  alter  breast  cancer  risk.  Understanding  these 
mechanisms  will  ultimately  facilitate  the  design  of  safe  and  effective  hormonal  chemoprevention 
regimens.  Moreover,  the  development  and  testing  of  such  regimens  will  be  greatly  facilitated  by 
the  identification  and  use  of  intermediate  molecular  endpoints  that  accurately  detect  changes  in 
the  breast  associated  with  changes  in  breast  cancer  risk.  We  have  chosen  to  exploit  the  defined 
relationship  between  parity  and  carcinogenesis  in  the  breast  to  generate  surrogate  endpoint 
biomarkers  for  changes  in  the  breast  associated  with  a  reduction  in  breast  cancer  risk.  While  our 
findings  do  not  address  whether  the  molecular  and  cellular  alterations  identified  in  this  study  are 
causally  related  to  parity-induced  protection  against  breast  cancer,  they  do  suggest  promising 
new  avenues  for  investigation.  We  believe  that  such  biomarkers  will  ultimately  prove  essential 
for  understanding  the  molecular  and  cellular  basis  of  parity-induced  protection  against  breast 
cancer,  and  for  the  rational  design  and  testing  of  hormonal  chemoprevention  regimens  aimed  at 
mimicking  this  naturally  occurring  protective  event. 
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EXPERIMENTAL  PROCEDURES 

Animals  and  Tissues 

FVB  mice  and  Sprague-Dawley  rats  were  housed  under  barrier  conditions  with  a  12  hr 
light/dark  cycle  and  access  to  food  and  water  ad  libitum.  Parous  rodents  were  generated  by 
mating  4- week-old  mice  or  9-week-old  rats.  Following  parturition,  animals  were  allowed  to 
lactate  for  21  days,  at  which  time  litters  were  weaned.  Parous  animals  underwent  28  days  of 
postlactational  involution  prior  to  sacrifice,  at  which  time  the  #3-5  mammary  glands  were 
harvested  and  snap  frozen,  as  were  those  of  age-matched  nulliparous  controls.  With  the 
exception  of  glands  used  for  whole  mounts,  the  lymph  nodes  within  gland  #4  were  removed. 

Whole  Mounts  and  Histology 

Number  4  mammary  glands  were  mounted  on  glass  slides,  fixed  overnight  in  neutral 
buffered  formalin,  and  transferred  to  70%  ethanol.  For  whole  mounts,  glands  were  rinsed  in 
water  for  5  min  and  stained  in  a  filtered  solution  of  0.2%  carmine  (Sigma)  and  0.5%  aluminum 
potassium  sulfate  for  1-3  days.  Glands  were  then  dehydrated  sequentially  through  70%,  90%, 
100%  ethanol  for  15  min.  each,  then  de-fatted  and  stored  in  methyl  salicylate.  For  histological 
analysis,  fixed  glands  were  blocked  in  paraffin,  sectioned,  and  stained  with  hematoxylin  and 
eosin. 

RNA  Isolation  and  Northern  Analysis 

Snap-frozen  tissue  was  homogenized  in  guanidine  thiocyanate  supplemented  with  7  |xl/ml 
2-mercaptoethanol,  and  RNA  isolated  by  centrifugation  through  cesium  chloride  as  previously 
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described  (80).  Equal  amounts  of  RNA  from  each  of  15-20  mice  or  10  rats  were  combined  for 
each  independent  pool.  Total  RNA  was  separated  on  a  1%  LE  agarose  gel,  and  passively 
transferred  to  Gene  Screen  (NEN).  Northern  hybridization  was  performed  per  manufacturer’s 
instructions  using  PerfectHyb  Plus  Hybridization  Buffer  (Sigma)  and  ^^P-labeled  cDNA  probes 
corresponding  to  Genbank  sequences  represented  on  the  Affymetrix  oligonucleotide  microarray 
Mu6500  Gene  Chip. 

Oligonucleotide  Microarray  Hybridization  and  Analysis 

Approximately  40  pg  of  total  pooled  RNA  from  each  sample  was  used  to  generate  cDNA 
and  biotinylated  cRNA  as  described  (81).  Hybridization  to  a  set  of  Affymetrix  Mu6500K 
microarrays  was  performed  per  manufacturer’s  instructions.  Following  washing  and  staining 
with  streptavidin-phycoerythrin,  chips  were  scanned  using  a  Hewlett-Packard  Gene  Array 
Scanner.  Grid  alignment  and  raw  data  generation  was  performed  using  Affymetrix  GeneChip 
3.1  software.  Raw  gene  expression  levels  were  scaled  and  normalized  data  sets  were  compared 
using  Affymetrix  algorithms  to  identify  differentially  expressed  genes. 

For  clustering  analysis,  probe  sets  for  differentially  expressed  genes  identified  in  the 
Mu6500  array  analysis  were  mapped  to  corresponding  probe  sets  on  Affymetrix  MGU74 
microarrays  via  the  Unigene  and  LocusLink  databases.  A  list  of  orthologous  probe  sets  on 
Affymetrix  RGU34  chip  was  generated  using  matches  obtained  by  querying  the  each 
differentially  expressed  gene  identified  in  the  mouse  against  the  Homologene  database 
(http://www.ncbi.nlm.nih.gov/HomoloGene/).  Of  the  resulting  matches,  only  curated  and 
calculated  reciprocal  best  match  Homologene  hits  were  selected  for  further  analysis.  Data  were 
scaled  such  that  the  mean  signal  intensity  was  equivalent  across  all  array  samples  excluding  the 
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top  and  bottom  2%  of  data  points.  Samples  were  standardized  prior  to  cluster  analysis  such  that 
the  median  expression  level  and  standard  deviation  for  each  gene  equaled  0  and  1 ,  respectively, 
across  the  set  of  samples  being  clustered.  The  data  was  then  filtered  to  include  only  probe  sets 
for  genes  shown  to  be  differentially  expressed  between  parous  and  nulliparous  mouse  samples 
(Table  1),  or  the  orthologues  for  these  genes  in  the  rat.  Cluster  software  was  used  to  generate 
hierarchical  clustering  trees,  which  were  visualized  using  Treeview  (M.  Risen; 
http://www.microarrays.org/software). 

In  Situ  Hybridization 

In  situ  hybridization  was  performed  as  described  (41).  Antisense  and  sense  probes  were 
sythesized  with  the  Promega  in  vitro  transcription  system  using  ^^S-UTP  and  ^^S-CTP  from  the 
T7  and  SP6  RNA  polymerase  promoters  of  a  PCR  template  containing  the  same  sequences  used 
for  Northern  hybridization  analysis. 
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FIGURE  LEGENDS 


Figure  1.  Parity-induced  morphological  changes  in  mice  and  rats.  Carmine-stained  whole 
mounts  (left  6  panels)  from  parous  involuted  from  Sprague-Dawley  rats  and  C57B1/6  and  FVB 
miee  each  display  increased  ductal  branching  as  compared  to  the  age-matehed  nulliparous  gland. 
Magnification  25X.  Histological  analysis  of  hematoxylin  and  eosin-stained  sections  (right  2 
panels)  demonstrates  that  the  proportion  of  epithelial  cells  is  similar  in  the  nulliparous  and 
parous  mammary  glands.  Magnification  400X. 

Figure  2.  Confirmation  of  parity-induced  changes  in  gene  expression.  Northern  hybridization 
analysis  of  gene  expression  for  candidate  genes  identified  by  microarray  analysis.  Expression 
levels  were  determined  for  independent  pools  of  mammary  gland  total  RNA,  each  of  which  was 
derived  from  15-20  age-matehed  nulliparous  (lanes  1-3)  or  parous  (lanes  4-6)  animals.  Genes 
downregulated  as  a  consequence  of  parity  ineluded  superoxide  dismutase  III  (SOD3),  leptin 
{Ob),  MUC18,  and  Mr,  a  homologue  of  the  Ig  Superfamily.  Genes  upregulated  by  parity 
included  carbonic  anhydrase  isozyme  II  ( CA2),  adipocyte  dijferentiation-related  protein 
(ADFP),  adenosine  deaminase  (ADA),  and  carboxyl  ester  lipase  (CEL).  Comparable  epithelial 
cell  content  in  nulliparous  and  parous  mammary  glands  is  demonstrated  by  equivalent  levels  of 
CKI8  expression  when  normalized  to  P-actin  or  28S  rRNA  loading  controls. 

Figure  3.  Differentially  expressed  genes  identified  by  microarray  analysis  reproducibly 
distinguish  between  nulliparous  and  parous  mammary  tissues  in  mice  and  rats.  (A)  Six 
independent  FVB  samples  (3  parous  and  3  nulliparous  pooled  samples;  8-10  mice  per  pool)  were 
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analyzed  on  Affymetrix  MGU74A  microarrays  and  clustered  based  on  the  expression  patterns  of 
genes  identified  as  being  differentially  expressed  in  an  independent  set  of  parous  and  nulliparous 
FVB  samples  analyzed  on  Affymetrix  Mu6500  microarrays.  (B)  Six  independent  Balb/c 
samples  (3  parous  and  3  nulliparous  pooled  samples;  3  mice  per  pool)  were  analyzed  on 
Affymetrix  MGU74A  microarrays  and  clustered  based  on  the  expression  patterns  of  genes 
identified  as  being  differentially  expressed  in  an  independent  set  of  parous  and  nulliparous  FVB 
samples  analyzed  on  Affymetrix  Mu6500  microarrays.  (C)  Six  Lewis  rat  samples  (3  parous  and 
3  nulliparous  pooled  samples;  3-4  animals  per  pool)  were  analyzed  on  Affymetrix  RGU34A 
microarrays.  Genes  identified  as  being  differentially  expressed  in  parous  and  nulliparous  mice 
based  on  the  analysis  of  Mu6500  microarray  data  were  mapped  via  Homologene  to  the  rat 
genome  in  order  to  perform  clustering  operations  as  above.  (D)  A  subset  of  differentially 
expressed  genes  associated  with  the  TGF-fi3  pathway  were  used  to  cluster  the  12  murine 
nulliparous  and  parous  samples  described  in  A  and  B.  (E)  A  subset  of  differentially  expressed 
genes  encoding  epithelial  growth  factors  were  used  to  cluster  the  6  rat  samples  described  in  C. 
Genes  listed  in  Table  1  that  were  not  represented  on  the  MGU74A  or  RGU34A  microarrays  were 
omitted  from  the  clustering  analysis. 

Figure  4.  Differential  expression  of  lymphoid  and  myeloid  markers  identify  parity-induced 
changes  in  cell  populations  within  the  mammary  gland.  (A)  Northern  analysis  of  three 
independent  pools  of  nulliparous  (lanes  1  -3)  or  parous  (lanes  4-6)  mammary  gland  total  RNA 
isolated  from  FVB  mice  demonstrates  parity-induced  markers  for  macrophages  (MPEGl),  B- 
lymphocytes  (KLC),  and  T-lymphocytes  (TDAG).  Hybridization  to  a  probe  for  fi-actin  is  shown 
as  a  loading  control.  (B)  Northern  hybridization  of  MME,  KLC  and  Eta-1  probes  to  total 
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mammary  gland  RNA  isolated  from  female  mice  at  the  indicated  developmental  stages. 
Expression  levels  are  compared  to  those  of  fi-actin  to  account  for  dilutional  effects  due  to  large- 
scale  increases  in  milk  protein  gene  expression  during  late  pregnancy  and  lactation.  (C)  In  situ 
hybridization  analysis  of  MME,  KLC  and  Eta-1  expression  in  the  mammary  gland  at  the 
indicated  developmental  stages.  Bright-field  (top)  and  dark-field  (bottom)  photomicrographs  of 
murine  mammary  gland  sections  hybridized  with  S-labeled  MME,  KLC  and  Eta-1  specific 
antisense  probes.  No  signal  over  background  was  detected  in  sections  hybridized  with  sense 
probes  (data  not  shown).  Exposure  times  were  identical  for  all  dark-field  photomicrographs  of 
each  gene  to  facilitate  comparison  of  gene  expression  levels.  Magnification  300x. 

Figure  5.  Parity-induced  downregulation  of  amphiregulin  and  pleiotrophin  expression  in  the 
mammary  gland.  (A)  Northern  hybridization  analysis  of  three  independent  pools  of  nulliparous 
(lanes  1-3)  or  parous  (lanes  4-6)  mammary  gland  total  RNA  isolated  from  FVB  mice 
demonstrates  decreased  expression  of  mitogenic  signaling  molecules,  fi-actin  is  shown  as  a 
loading  control.  (B)  Northern  hybridization  analysis  of  the  developmental  patterns  of  Areg  and 
Ptn  expression  in  the  mammary  gland.  Expression  levels  are  compared  to  those  of  |3-actin  to 
account  for  dilutional  effects  due  to  large-scale  increases  in  milk  protein  gene  expression  during 
late  pregnancy  and  lactation.  The  28S  rRNA  band  is  shown  as  a  loading  control.  (C)  In  situ 
hybridization  analysis  of  Areg  expression  at  the  indicated  developmental  stages.  Bright-field 
(top)  and  dark-field  (bottom)  photomicrographs  of  murine  mammary  gland  sections  hybridized 
with  an  ^^S-labled  ^rcg-specific  antisense  probe.  No  signal  over  background  was  detected  in 
sections  hybridized  with  a  sense  Areg  probe  (data  not  shown).  Exposure  times  were  identical  for 
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all  dark-field  photomicrographs  to  facilitate  comparison  of  gene  expression  changes. 
Magnification  300x. 

Figure  6.  Parity-induced  increases  in  TGF-fiS  and  clusterin  expression  in  the  mammary  gland. 
(A)  Northern  analysis  of  three  independent  pools  of  nulliparous  (lanes  1-3)  and  parous  (lanes  4- 
6)  mammary  gland  RNA  isolated  from  FVB  mice  demonstrates  increased  expression  of  TGF-^3 
and  clusterin  in  the  parous  gland,  fi-actin  is  shown  as  a  loading  control.  (B)  Northern 
hybridization  analysis  of  the  developmental  patterns  of  TGF-fiS  and  clusterin  expression. 
Hybridization  to  ji-actin  is  shown  s  control  for  loading  and  for  dilutional  effects  of  milk  protein 
gene  expression.  (C)  In  situ  hybridization  analysis  of  TGF-fiS  and  clusterin  expression  at  the 
developmental  stages.  Bright-field  (top)  and  dark-field  (bottom)  photomicrographs  of  murine 
mammary  gland  sections  hybridized  with  ^^S-labeled  TGF-fiS  or  clusterin  antisense  probe. 
Magnification  300x. 

Figure  7.  Parity-induced  changes  in  gene  expression  are  permanent  and  conserved  across 
species.  A)  Northern  hybridization  analysis  for  the  indicated  genes  performed  on  pools  of  RNA 
from  the  mammary  glands  of  parous  animals  that  had  undergone  increasing  periods  of  involution 
and  of  their  age-matched  nulliparous  controls.  Parous  animals  were  mated  at  4  wks  of  age,  and 
underwent  21  days  of  lactation  and  either  4  wks,  16  wks,  or  30  wks  of  postlactational  involution. 
B)  Probes  for  the  indicated  genes  were  hybridized  to  Northern  membranes  containing  pools  of 
total  RNA  from  3  parous  (P)  and  3  age-matched  nulliparous  (N)  129SvEv  and  Balb/c  mice  C) 
Probes  for  the  indicated  rat  genes  were  hybridized  to  Northern  membranes  containing  pools  of 
total  RNA  from  10  parous  (P)  and  10  nulliparous  (N)  Sprague-Dawley  rats.  In  each  case,  both 
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the  direction  and  magnitude  of  change  in  expression  between  the  nulliparous  and  parous  states 
are  comparable  to  that  observed  in  FVB  mice.  P-actin  and  28S  rRNA  are  shown  as  a  loading 
control. 
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Table  1:  Differentially  Expressed  Genes 


Genes  Down-Regulated 


Sequence  Identity 

Accession  Number 

Function 

Avg.  Fold  Change 

Northern 

Amphiregulin 

L41352 

Growth  Factor 

2.8 

Y 

Pleiotrophin 

D90225 

Growth  Factor 

2.0 

Y 

Insulin-Like  Growth  Factor  1B 

W10072 

Growth  Factor 

2.2 

ND 

Insulin-Like  Growth  Factor  1A 

X04480 

Growth  Factor 

1.6 

Y 

Thyroid  Stimulating  Hormone  Receptor 

U02602 

Hormone  Receptor 

2.0 

Y 

Leptin 

U18ai2 

Hormone 

2.0 

Y 

Ig  Superfamily  Containing  Leucine-rich  repeat 

AA059664 

Cell  Adhesion 

2.5 

Y 

MUC  18 

AA088962 

Ceil  Adhesion 

2.0 

Y 

Superoxide  Dismutase  III  (SOD3) 

X84940 

Oxidoreductase 

2.5 

Y 

Carbonic  Anhydrase  III 

M27796 

Hydratase 

1.4 

Y 

Genes  Up-Regulated 


Sequence  Identity 

Accession  Number 

Function 

Avg.  Fold  Change 

Northern 

Whey  Acidic  Protein 

J00544 

DifferentiationA4itk  Protein 

24 

Y 

Gamma-Casein 

D10215 

Oifferentiatlon/Miik  Protein 

20 

Y 

Alpha-Lactalbumin 

M80909 

Differentiation/Mllk  Protein 

13 

Y 

Lactoferrin 

J03298 

Differenliation/lron  Transport 

4.8 

Y 

Alpha-Casein 

M36780 

Differenliation/Milk  Protein 

4.2 

ND 

WDNM1 

X93037 

Differentiation/Protease  Inhibitor 

3.5 

Y 

Beta-Casein 

X04490 

Differentiation/Milk  Protein 

3.3 

Y 

Kappa-Casein 

M10114 

Oifferentiation/Milk  Protein 

2.9 

Y 

Adipocyte  Differentiation  Related  Protein 

M93275 

Differentiation 

1.8 

Y 

Carboxyl  Ester  Lipase 

U37386 

Lipid  Degradation 

6.6 

Y 

LPS-Binding  Protein 

X99347 

Antibacterial/Milk 

6.5 

Y 

Lysozyme  P 

M21050 

Antibacterial/Milk 

2.0 

ND 

Immunoglobulin  M  Heavy  Chain 

ET61785 

Immunoglobulin 

30 

Y 

Immunoglobulin  G  Heavy  Chain 

ET61798 

Immunoglobulin 

13 

Y 

Immunoglobulin  A  Heavy  Chain 

J00475 

Immunoglobulin 

21 

Y 

Kappa  Light  Chain 

XI 6678 

Immunoglobulin 

7.0 

Y 

Macrophage  Metalloelastase 

M82831 

Metalloprotease 

2.1 

Y 

Macrophage  Expressed  gene  1  (Mpeg1) 

L20315 

Cell  Signalling 

1.7 

Y 

Lipocalin  2 

W13166 

Cell  Signalling 

2.6 

Y 

T-cell  death  associated  gene  (Tdag) 

U44088 

Celt  Signalling 

2.3 

Y 

Transforming  Growth  Factor  Beta-3 

M32745 

Growth  Inhibition 

1.9 

Y 

Clusterin 

L08325 

Apoptosis/TGF-(3  Pathway 

2.0 

Y 

Eta-1 

X16151 

Cell  Signalling/TGF-B  Pathway 

8.8 

Y 

Id-2 

M69293 

Cell  Cycle/TGF-B  Pathway 

2.3 

Y 

Carbonic  Anhydrase  II 

K00811 

Hydratase 

3.0 

Y 

Cyclin  D1 

A1849928 

Cell  Cycle 

1.8 

Y 

CRBP1 

X60367 

Retinol-Binding 

3.5 

Y 

Connexin  26 

M81445 

Gap  Junction 

3.1 

Y 

Folate-Binding  Protein  1 

ET63126 

Receptor 

2.6 

Y 

Adenosine  Deaminase 

M10319 

Nucleotide  Metabolism 

3.0 

Y 

Chitinase  3-like  1  (BRP-39) 

X93035 

Glycoprotein 

3.2 

Y 

Table  1:  Genes  found  to  exhibit  a  persistent  change  in  expression  in  fully  involuted  compared  to  age-matched 
nulliparous  mammary  glands.  Genes  whose  differential  expression  was  confirmed  by  Northern  analysis  are 
indicated  (Y).  Genes  whose  differential  expression  was  not  tested  (ND)  by  Northern  hybridization  were  listed 
if  they  displayed  a  consistent  change  in  regulation  in  3/3  microarray  experiments. 
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MEMORANDUM  FOR  Administrator,  Defense  Technical  Information 
Center  (DTIC-OCA) ,  8725  John  J.  Kingman  Road,  Fort  Belvoir, 

VA  22060-6218 

SUBJECT;  Request  Change  in  Distribution  Statement 

1.  The  U.S.  Army  Medical  Research  and  Materiel  Command  has 
reexamined  the  need  for  the  limitation  assigned  to  technical 
reports  written  for  this  Command.  Request  the  limited 
distribution  statement  for  the  enclosed  accession  numbers  be 
changed  to  "Approved  for  public  release;  distribution  unlimited 
These  reports  should  be  released  to  the  National  Technical 
Information  Service. 


2.  Point  of  contact  for  this  request  is  Ms.  Kristin  Morrow  at 
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